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Abstract 

Decontamination  activities  may  cause  the  release  of  contaminated  washwater  into 
the  wastewater  that  eventually  flows  into  a  wastewater  treatment  facility.  This  raises 
concerns  about  the  effect  of  chemical  warfare  agents  (CWAs)  on  the  microbial  con¬ 
sortia  that  are  responsible  for  cleaning  the  wastewater.  This  study  investigated  the 
impact  of  malathion  on  the  microbial  ecology  of  laboratory  scale  activated  sludge  com¬ 
munities.  The  Simpson  Reciprocal  Index  decreased  for  three  bioreactors  operated  in 
the  absence  of  malathion,  which  showed  that  the  microbial  assemblage  became  less 
diverse  during  the  course  of  the  study.  The  species  identified  in  the  bioreactors  be¬ 
longed  to  well-known  groups  of  heterotrophic  and  autotrophic  bacteria,  and  these 
groups  were  represented  in  the  presence  and  absence  of  malathion.  Nitrospira,  a 
key  player  in  autotrophic  nitrogen  removal,  decreased  in  relative  abundance  for  the 
bioreactor  exposed  to  0.1  ^  of  malathion  but  increased  in  the  bioreactor  exposed 
to  3  ^  of  malathion,  possibly  due  to  interactions  with  heterotrophic  groups  that 
suffered  inhibition.  To  the  author’s  knowledge,  this  study  is  the  Erst  to  document  the 
ecological  impacts  of  long  term  malathion  exposure  in  bioreactors  carrying  out  COD 
removal  and  nitrification. 
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EFFECT  OF  MALATHION  ON  THE  MICROBIAL  ECOLOGY  OF  ACTIVATED 


SLUDGE 

I.  Introduction 

1.1  Motivation  and  Background 

Organophosphates  (OPs)  are  widely  used  in  the  application  of  pesticides.  Fur¬ 
thermore,  some  countries  persist  in  maintaining  chemical  warfare  stocks  as  insurance 
against  neighboring  states  with  similar  capabilities  or  more  technologically  advanced 
conventional  weapons.  Non-state  actors  continue  to  explore  weapons  that  generate 
publicity  and  differentiate  their  organization  from  others  with  more  conventional  ca¬ 
pabilities.  The  public  risk  associated  with  exposure  to  OPs  demands  an  adequate 
response  to  mitigate  the  hazard.  This  typically  would  involve  a  large  quantity  of  wa¬ 
ter  and  a  strong  alkaline  or  chlorinated  neutralizing  agent,  which  would  be  collected 
and  treated  as  hazardous  waste  or  in  a  wastewater  treatment  facility. 

Most  municipal  wastewater  treatment  plants  in  the  United  States  use  activated 
sludge  systems.  Activated  sludge  systems  utilize  a  process  that  includes  filtration, 
sedimentation  and  biodegradation  to  eliminate  contaminants  from  an  effluent  stream. 
Early  phases  remove  large  debris  through  filtration.  Secondary  phases  allow  larger 
particulates  to  settle  in  a  sedimentation  tank.  Bioremediation  occurs  in  an  aerobic 
reactor  where  activated  sludge  containing  nitrifying  bacteria  and  other  microorgan¬ 
isms  break  down  pollutants  by  feeding  on  oxidation  products  from  compounds  in 
the  sludge.  A  final  treatment  stage  removes  contaminants  using  filtration  before  the 
treated  effluent  is  released  into  the  environment  (1). 
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A  diverse  community  of  microbes  develops  using  a  suitable  feedstock  in  the  biodegra¬ 
dation  phase  to  reduce  or  eliminate  the  contaminant.  An  aerobic  sequencing  batch 
reactor  is  an  environment  where  the  nitrogen  cycle,  a  crucial  process  in  the  biosphere, 
is  utilized  to  transform  nitrogen  compounds  into  energy  sources  for  bacteria.  Wastew¬ 
ater  treatment  plants  take  advantage  of  the  nitrogen  cycle  and  other  transformation 
processes  to  biodegrade  toxic  substances  into  less  harmful  substances. 

Recent  studies  indicate  that  the  composition  of  the  community  may  contribute  to 
the  effectiveness  of  its  biodegradation  capability  with  respect  to  OPs  (2)  (3).  Fur¬ 
ther  study  is  necessary  to  identify  which  genera  are  most  able  to  survive  in  an  OP 
environment  and  whether  specific  genera  or  combinations  may  be  more  effective  at 
biodegradation.  Various  methods  exist  for  the  evaluation  of  bacterial  communities.  In 
recent  years,  the  development  of  next  generation  sequencing  technology  has  extended 
the  breadth  and  depth  of  evaluation  results.  The  application  of  sequencing  technology 
offers  further  phylogenetic  detail  that  may  improve  bioremediation  research  efforts. 
Bioremediation  offers  an  economical  and  efficient  means  to  treat  water  supplies  and 
developing  optimal  parameters  ensures  effective  elimination  of  OP  contaminants. 

1.2  Research  Objective 

This  research  project  inquires:  What  is  the  impact  of  malathion  to  bacterial  con¬ 
sortia  in  activated  sludge  systems?  Which  bacteria  are  most  impacted  by  the  presence 
of  malathion?  What  traits  are  displayed  by  bacteria  that  exist  in  environments  con¬ 
taminated  with  malathion?  The  objectives  for  the  research  are  assess  the  diversity  of 
the  bacterial  community  and  determine  any  effects  on  the  community  structure. 

In  order  to  simulate  the  environment  where  an  OP  may  be  introduced  into  a 
wastewater  treatment  plant,  it  was  necessary  to  build  three  bench-scale  activated 
sludge  sequencing  batch  reactors.  The  setup  consisted  of  three  containment  vessels 


2 


with  2.0  L  capacity,  each  with  different  concentrations  of  substrate.  The  three  reactors 
were  seeded  with  activated  sludge  from  a  local  wastewater  treatment  plant  and  the 
sludge  was  provided  with  simulated  wastewater  containing  macro  and  micro  nutrients. 
Malathion,  a  commonly  used,  commercially  available,  and  less  toxic  OP  pesticide  was 
used  as  a  surrogate  for  higher  toxicity  nerve  agents.  After  a  period  of  acclimation  to 
the  laboratory  environment,  malathion  was  added  to  the  influent  and  pumped  into  two 
of  the  reactors  continuously  for  30  days  maintaining  reactor  substrate  concentrations 
of  0.1  Y  and  3  rf. 

The  study  evaluated  the  ecology  of  bacteria  when  exposed  to  OPs.  Genomic  DNA 
from  biomass  samples  was  extracted  using  procedures  adopted  from  Janeczko  (1). 
The  MoBio  PowerSoil  DNA  isolation  kit  was  used  to  extract  nucleic  material  from 
bacterial  samples  for  PCR  analysis  and  genetic  identification.  The  three  reactors 
were  evaluated  using  genetic  sequencing  over  an  extended  period  prior  to  malathion 
exposure  to  establish  a  baseline  phylogeny  and  identify  community  changes  due  to  the 
transition  to  a  laboratory  environment.  Additional  samples  were  taken  after  exposure 
for  a  period  of  one  month  and  sequenced  using  next  generation  sequencing. 
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II.  Literature  Review 


2.1  Organophosphates 

Organophosphates  have  long  been  used  in  applications  as  pesticides  and  military 
nerve  agents.  The  toxicity  of  OPs  is  well-known  and  contributed  to  their  historical 
development  as  weapons  of  mass  destruction  in  the  inventories  of  Nazi  Germany,  So¬ 
viet  Russia  and  the  United  States  (4).  Common  examples  of  pesticides  are  malathion, 
parathion,  and  diazinon.  Perhaps  more  well-known  are  the  notorious  chemical  nerve 
agents  VX,  sarin,  soman,  and  tabun.  The  toxicity  of  OPs  varies  with  exposure  and 
the  species  of  OP  encountered.  Poisoning  cases  associated  with  OPs  occur  worldwide 
and  are  frequently  associated  with  suicides  in  agrarian  locations  (5).  An  epidemic 
poisoning  occurred  in  Pakistan  where  7,500  workers  were  exposed  to  malathion,  an 
OP  insecticide.  At  its  peak,  it  is  estimated  that  2,800  workers  were  casualties  and 
at  least  five  died  from  the  effects  (6).  Malathion  degrades  to  more  toxic  products 
such  as  malaoxon  and  isomalathion.  Malaoxon  is  an  oxidative  metabolite  that  is 
61  times  more  toxic  than  malathion.  Isomalathion  is  1000  times  more  toxic  (7). 
However,  mammals  and  birds  have  higher  carboxylesterase  activity  than  pests  such 
as  insects,  allowing  malathion  removal  at  lower  concentrations  and  preventing  more 
toxic  metabolites  from  accumulating  (8). 

Early  research  into  OP  chemistry  and  use  began  in  the  1800s  with  French  organic 
chemist  P.  Clermont  synthesizing  the  first  cholinesterase  inhibitor,  TEPP.  In  1932, 
a  German  chemist  first  described  nervous  systems  effects  including  a  choking  sensa¬ 
tion  and  dimming  of  vision.  Another  German  chemist,  G.  Schrader,  studied  OPs  in 
hopes  of  producing  new  insecticides,  but  later  became  known  as  “the  father  of  nerve 
agents”  for  his  accidental  discovery  of  both  sarin  and  tabun.  The  Nazi  government 
seized  on  his  discovery,  placing  him  in  charge  of  a  chemical  warfare  agent  program 
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that  produced  large  quantities  of  OPs  called  G-agents.  Another  series  of  chemical 
agents,  the  V-agents,  was  first  produced  by  the  British  during  World  War  II.  After 
the  war,  American  companies  used  Schrader’s  research  to  produce  pesticides  includ¬ 
ing  malathion,  parathion  and  azinphosmethyl.  OP  use  as  a  pesticide  was  boosted 
later  when  organo chlorine  pesticides  were  banned  in  the  1970s  (9). 

Organophosphates  are  esters  of  phosphoric  acid  with  applications  as  pesticides 
and  are  characterized  by  their  chemical  structural  features.  These  include  a  terminal 
oxygen  connected  to  a  phosphorus  by  a  double  bond.  The  phosphorus  is  bonded  to 
two  lipophilic  groups  and  a  leaving  group,  which  is  often  a  halide.  Examples  of  this 
type  of  OP  are  the  nerve  agents  VX  and  sarin.  Some  compounds  replace  the  double- 
bonded  oxygen  with  a  terminal  sulfur.  Pesticides  such  as  parathion  and  malathion 
have  this  feature.  In  mammals,  metabolism  removes  a  lipophilic  group,  resulting  in 
toxic  effects  less  severe  than  that  of  nerve  agents.  In  insects,  metabolism  involves 
oxidation  of  the  compound,  replacing  the  sulfur  with  a  terminal  oxygen  improving  its 
enzyme  inhibition  capability  (9). 

However,  both  groups  function  in  a  similar  manner.  In  the  human  body,  OPs  enter 
through  dermal,  inhalation,  or  ingestion  routes  and  are  distributed  throughout  the 
body.  When  OPs  reach  the  nervous  system,  they  have  significant  toxic  effects.  OPs 
inhibit  acetylcholinesterase  (AChE),  an  enzyme.  Inhibition  is  caused  by  phosphoryla¬ 
tion  of  the  serine  hydroxyl  residue,  inactivating  AChE.  This  leads  to  the  accumulation 
of  acetylcholine  (ACh)  at  postsynaptic  membranes.  In  early  stages,  the  phosphoryla¬ 
tion  may  be  reversible,  but  eventually  aging  occurs  resulting  in  permanent  covalent 
bonds.  Excess  ACh  causes  permanent  stimulation  at  cholinergic  nerve  junctions  and 
in  the  muscles.  Symptoms  may  progress  from  parasympathetic  and  sympathetic  hy¬ 
peractivity  to  lethal  effects  when  OPs  are  present  in  sufficient  concentration  (9). 
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2.2  Degradation 


Numerous  studies  have  documented  various  means  for  the  biodegradation  of  pes¬ 
ticides  and  chemical  warfare  agents.  Insecticides  generally  fall  within  one  of  three 
categories  of  organophosphates:  phosphotriesters,  thiophosphatetriesters,  or  phos- 
phorothioesters.  The  primary  difference  between  each  is  the  number  of  oxygen  or 
sulfur  groups  attached  to  the  phosphorus  center.  In  the  case  of  malathion,  a  phos- 
phorothioester,  the  phosphorus  center  is  double  bonded  to  a  sulfur,  with  an  opposing 
sulfur  leaving  group  bonded  to  the  phosphorus  center  (10).  These  varying  struc¬ 
tures  play  a  role  in  the  effectiveness  of  microorganisms’  mechanisms  to  utilize  the 
constituent  chemical  bonds  to  satisfy  energy  and  nutrition  requirements. 

The  primary  tool  that  bacteria  and  other  microorganisms  use  is  enzymes.  Enzy¬ 
matic  hydrolysis  generally  occurs  using  a  multistep  processes  involving  phosphatase, 
esterase,  hydrolase  and  oxygenase  (11).  Some  enzymes  have  broad  substrate  speci¬ 
ficity,  such  as  organophosphate  hydrolase  (OPH),  which  is  used  by  bacteria  like  Pseu¬ 
domonas  diminuta  (12).  OPH  has  shown  the  capability  to  degrade  a  wide  variety  of 
OPs  into  less  toxic  end  products.  It  has  also  been  shown  to  be  an  effective  enzyme  in 
degrading  the  highly  toxic  VX  nerve  agent,  along  with  other  phosphate-sulfur  bond 
OPs  such  as  demoton,  malathion,  acephate,  and  phosalone  (13)  (14).  Often,  one  or 
more  enzymes  is  required  and  may  produce  varying  levels  of  results  depending  on 
the  substrate  involved.  Another  widely  studied  enzyme  is  organophosphorus  acid 
anhydrolase  (OPAA),  isolated  from  Alteromonas.  OPAA  works  well  at  hydrolyzing 
phosphorus-flourine  bonds,  but  works  slowly  against  phosphorus  bonds  with  oxygen 
and  doesn’t  work  against  sulfur  (15). 

Several  species  of  soil  bacteria  including  Bacillus,  Pseudomonas,  Flavobacterium , 
Arthrobacter ,  Micrococcus  and  others  have  shown  varying  capability  to  degrade  OPs 
through  species-unique  enzymatic  processes  (16).  Bacillus  shows  significant  phos- 


6 


phatase  activity,  particularly  at  higher  pH  levels  (~8.0)  (17).  Pseudomonas  sp.  are 
frequently  named  as  a  successful  genera  capable  of  degrading  malathion  using  en¬ 
zymatic  hydrolysis  via  dicarboxylesterase  or  carboxylesterase  (18).  Degradation  by 
Rhizobium  sp.  has  also  been  reported  implicating  carboxylesterase  as  the  likely  path¬ 
way  (19).  Some  organisms  require  additional  carbon  sources  in  cometabolism  pro¬ 
cesses  to  successfully  degrade  organophosphorous  compounds.  Sodium  succinate  and 
sodium  acetate  are  examples  of  cometabolism  compounds  that  enhance  malathion 
degradation  by  Acinetobacter  johnsonii  MA19  (20).  Furthermore,  a  genetic  basis  for 
OP  degradation  has  been  established.  Organophosphorous  degrading  bacteria  have 
been  genetically  engineered  with  genes  that  express  enzymatic  capability  and  proven 
successful.  In  one  example  of  genetically  engineered  bacteria,  Pseudomona  putida 
was  isolated  from  p-nitrophenol  contaminated  waste  and  engineered  to  degrade  OP 
compounds  simultaneously  with  p-nitrophenol  (21).  Additionally,  there  is  some  evi¬ 
dence  to  suggest  plasmid  involvement  during  malathion  degradation  in  Micrococcus 
sp.  (12). 

2.3  Waste  Processing 

Wastewater  treatment  involves  a  series  of  filtering,  sedimentation,  and  treatment 
processes  that  prevent  undesirable  substances  from  reemerging  in  water  supplies  for 
human  consumption.  Initially,  wastewater  passes  through  screens  that  remove  large 
debris.  It  then  flows  to  a  settling  tank  where  solid  particles  are  allowed  to  precipitate 
and  the  fluid  portion  is  separated  from  it.  The  fluid  is  aerated  and  a  nitrifying 
bacterial  community  develops  in  the  nutrient-rich,  oxidation  tank.  At  this  point,  the 
mixture  is  called  activated  sludge. 

Activated  sludge  consists  of  a  wide  variety  of  microorganisms.  One  group  of 
bacteria  involved  in  the  biodegradation  process  is  the  nitrifying  bacteria.  Nitrify- 
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ing  bacteria  oxidize  inorganic  nitrogen  compounds  for  energy  consumption.  These 
can  be  further  subdivided  into  ammonia  oxidizers  and  nitrite  oxidizers.  Ammonia 
oxidizers,  including  members  from  Proteobaeteria,  are  represented  by  genera  Nitro- 
somonas,  Nitrosospira,  Nitrosolobus,  and  Nitrosovibrio  (beta-class)  and  Nitrosococcus 
(gamma-class).  The  nitrite  oxidizers  are  also  members  of  Proteobaeteria,  including 
genera  Nitrobacter  (alpha-class),  Nitrococcus  and  Nitrospina  (both  gamma-class). 
Nitrospira  ( Nitrospira  phylum)  is  also  a  nitrite  oxidizer.  The  ammonia  oxidizers  are 
lithotrophic  bacteria  that  oxidize  ammonia  into  nitrite  through  the  enzyme  ammonia 
monooxygenase  (AMO),  and  then  the  nitrite  oxidizers  oxidize  the  nitrite  into  nitrate 
using  the  enzyme  nitrite-oxidoreductase  (NO2-OR)  (17). 

The  microorganisms  are  fed  a  mixture  of  nutrients  and  begin  consuming  dissolved 
compounds  in  the  solution.  After  a  certain  period  of  time,  the  sludge  settles  to 
the  bottom  where  it  may  be  recycled  for  use  again,  disposed  in  effluent  wastewater 
streams,  or  sent  to  a  waste  processing  area  for  disposal  at  a  landfill  or  as  fertilizer. 
The  amount  of  time  spent  in  the  activated  sludge  phase  affects  the  removal  rate  of 
contaminants.  The  contaminant  removal  capacity  of  the  activated  sludge  is  depen¬ 
dent  on  factors  such  as  the  microorganism  community  structure  and  performance. 
Additionally,  since  the  activated  sludge  is  a  biomass,  the  contaminant  profile  and 
concentration  of  its  constituents  will  affect  the  ability  of  living  organisms  exposed 
to  it  to  perform  their  intended  functions  (1).  Understanding  the  effects  of  OP  con¬ 
taminants  and  the  ability  of  a  wastewater  treatment  facility  to  effectively  eliminate 
their  toxicity  is  a  key  component  of  the  final  phase  of  recovery.  Once  waste  leaves  the 
facility,  it  will  be  reintroduced  to  the  environment  where  it  may  once  again  present 
potential  exposure  risk. 


2.4  Phylogenetic  Analysis  of  Microbial  Communities 


Microbial  detection  and  community  structure  phylogeny  are  determined  using  var¬ 
ious  methods.  Bergey’s  Manual  of  Determinative  Bacteriology  provides  a  flow  chart 
method  for  identification  and  comparison  of  observed  sample  culture  characteristics 
with  known  bacterial  characteristics.  Bergey’s  Manual  was  used  to  identify  the  Pseu¬ 
domonas  genus  during  an  evaluation  of  malathion  degradation  during  continuous 
cultivation  in  diluted  and  undiluted  nutrient  broths  (16).  This  analytical  approach 
involves  morphological  and  cultural  means  to  identify  the  organism’s  physiology. 

In  the  last  few  decades,  molecular  techniques  that  amplify  nucleic  material  us¬ 
ing  polymerase  chain  reaction  (PCR)  have  become  more  common  and  specific.  The 
replication  process  for  PCR  involves  heating  and  denaturing  the  DNA  to  unwind  the 
double  helix  and  separate  the  two  sugar-phosphate  backbones  by  disassociating  the 
hydrogen  bonds  between  the  nucleotide  base  pairs.  The  reaction  then  cools  during 
the  annealing  step,  which  allows  a  forward  and  reverse  primer  to  attach  and  bind 
to  the  target  region  at  primer-specific  temperatures.  The  reaction  continues  with  a 
rise  to  the  optimal  polymerase  temperature  so  that  the  enzyme  can  synthesize  and 
elongate  the  new  DNA  strand  with  deoxynucleotide  triphosphate  (dNTPs)  bases  that 
have  been  added  to  the  reaction.  The  goal  of  PCR  is  to  produce  exponentially  greater 
quantities  of  the  target  region,  so  the  cycle  is  repeated  several  times. 

The  PCR  process  was  further  improved  and  evaluated  using  DNA  sequencing. 
DNA  sequencing  is  the  process  of  analyzing  DNA  strands  to  determine  the  order  of 
the  constituent  nucleotides.  In  this  study,  a  relatively  new  technology  called  massively 
parallel,  or  next  generation  sequencing  (NGS),  is  used  to  evaluate  samples  for  bacterial 
community  structure.  Massively  parallel  sequencing  is  possibly  more  descriptive, 
alluding  to  the  computing  power  required  to  identify  nucleotide  bases  real-time  during 
PCR.  All  NGS  instruments  include:  random  fragmentation  and  ligation  with  a  custom 
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library,  amplification  on  a  bead  or  glass,  direct  detection  during  sequencing,  millions 
of  reactions  during  a  run,  shorter  reads  than  capillary  sequencers,  and  digital,  paired- 
end  (from  both  ends)  reads  for  more  accurate,  mapped  quantitative  analysis  (22). 

The  16S  ribosomal  ribonucleic  and  deoxyribonucleic  acid  gene  sequencing  analysis 
method  is  one  of  the  more  common  methods  of  analyzing  bacterial  phylogeny.  The 
16S  gene  is  present  in  almost  all  bacteria  and  has  been  conserved  over  time  (23).  It 
is  also  specific  enough  for  even  rare  isolates  due  to  its  ubiquitous  presence,  size  and 
functional  purpose  (24).  Alternative  quantitative  techniques  are  also  available,  such 
as  automated  ribosomal  intergenic  spacer  analysis  (ARISA),  which  provided  a  com¬ 
munity  profile  as  an  exploratory  tool  in  the  case  of  Candidatus  Accumulibacter  (25). 
The  16S  method  is  used  in  this  study. 

2.5  Conclusion 

In  summary,  organophosphates  have  existed  for  many  years,  but  only  in  the  last 
two  centuries  have  their  chemical  properties  and  toxic  effects  become  well-known.  The 
study  of  OPs  has  led  to  the  creation  of  military  nerve  agents  and  pesticides  that  are 
still  in  use.  Organophosphates  bind  to  ACh  enzymes  and  cause  overstimulation  at  the 
post-synaptic  junctions  between  the  nervous  system  and  muscles.  Overstimulation 
may  be  fatal.  When  OP  hazardous  materials  incidents  or  terrorist  attacks  occur,  the 
response  will  likely  involve  decontamination  and  disposal.  Disposal  at  wastewater 
treatment  facilities  is  an  economical  and  possibly  unavoidable  consequence. 

Most  municipal  wastewater  treatment  plants  in  the  United  States  use  activated 
sludge  systems.  Activated  sludge  systems  utilize  a  process  that  includes  filtration, 
sedimentation  and  biodegradation  to  eliminate  contaminants  from  an  effluent  stream. 
A  diverse  community  of  microbes  develops  using  a  suitable  feedstock  in  the  biodegra¬ 
dation  phase  to  reduce  or  eliminate  the  contaminant.  Recent  studies  indicate  that  the 
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composition  of  the  community  may  contribute  to  the  effectiveness  of  its  biodegrada¬ 
tion  capability  with  respect  to  OPs.  This  community  may  be  characterized  in  greater 
detail  using  the  recently  developed  sequencing  technologies. 
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III.  Methodology 


3.1  Laboratory  Preparation 

An  initial  assessment  of  laboratory  conditions  and  resources  was  conducted  one 
month  prior  to  obtaining  the  sample  from  the  water  treatment  plant.  An  initial 
inventory  of  equipment  included  a  2  L  bench-scale  sequencing  batch  reactor,  peri¬ 
staltic  pumps  and  associated  hoses  and  tubing,  gas  and  liquid  solenoids,  electronically- 
controlled  switches,  laptops  with  control  software,  liquid  reservoir  and  containment 
trays,  and  associated  mounting  hardware.  Each  equipment  item  was  inspected  for 
cleanliness  and  operational  condition.  A  cleaning  solution  consisting  of  Alconox  and 
deionized  water  was  used  to  clean  the  tubing  and  reactors  and  then  rinsed  with  deion¬ 
ized  water  until  all  remaining  residue  was  removed.  All  reactor  pumps  and  lengths  of 
tubing  were  calibrated  using  15  runs  of  deionized  water  to  verify  volumetric  through¬ 
put. 

The  three  continuously  operated  benchtop  reactors  were  assembled  (Fig.  1).  An 
influent  water  tank  and  two  feed  solutions  (described  below)  in  1  L  (A)  and  0.5  L 
(B)  glass  containers  with  plastic  screw-on  lids  were  connected  via  tubing  through  a 
peristaltic  pump  to  each  of  the  three  reactors.  A  glass  tube  extending  from  each  of 
the  reactor  lids  to  the  liquid  surface  was  connected  via  tubing  through  an  electric 
solenoid  valve  to  the  effluent  tank.  Compressed  air  at  1-3  psi  was  pumped  into  each 
tank  via  tubing  through  an  electric  solenoid  valve.  Filtered  air  in  and  air  out  were 
also  connected  to  each  tank  via  tubing  and  controlled  with  an  electric  solenoid  valve. 
Control  of  all  pumps  and  switches  was  maintained  using  a  laptop  and  a  software- 
controlled  timing  switch.  Each  of  the  reactors  were  isolated  from  the  others  with  no 
transference  of  air  or  liquid  between  them. 

A  feed  solution  was  prepared  using  methods  adopted  from  previous  studies  in 
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Figure  1.  Reactor  Diagram 

the  lab  (1).  Stock  solutions  were  prepared  as  follows.  Feed  A  consisted  of  89.2  g 
NaHC03  and  2  L  of  deionized  water.  Feed  B  consisted  of  2.5  g  C2H3Na02,  12.0  g 
Peptone,  4.52  g  NH4C1,  13.72  g  MgCl2«6H20,  3.44  g  CaCl2«2H20,  1.335  g  KH2P04, 
40  mL  trace  element  solution  and  2  L  deionized  water.  The  trace  element  solu¬ 
tion  consisted  of  2.73  g  C6H807,  2.00  g  C9H9N03,  0.36  g  C6H6NNa306»2H20,  0.15 
g  CioHi3N2Na3Os»4H20,  1.50  g  FeCl3«6H20,  0.25  g  H3B03,  0.15  g  ZnS04«7H20, 
0.12  g  MnCl2«4H20,  0.07  g  CuS04»5H20,  0.03  g  KI,  0.03  g  Na2Mo04«2H20,  0.03  g 
CoC12«6H20,  0.03  g  NiCl2«6H20,  0.03  g  Na2W04«2H20,  and  1  L  of  deionized  water. 
All  three  stock  solutions  were  stored  in  a  laboratory  refrigerator  and  dispensed  as 
needed  into  autoclaved  glass  containers  for  use. 

3.2  Sample  Seed 

The  activated  sludge  was  obtained  from  Fairborn,  Ohio  Water  Treatment  Facility 
on  Day  1.  A  4  L  sample  was  removed  from  an  oxidation  tank  and  funncled  into  a 
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glass  screw-top  container.  Heavy  seasonal  precipitation  contributed  to  higher  than 
normal  dilution  of  the  sample.  The  4  L  sample  was  transported  to  the  laboratory 
and  allowed  to  settle.  2  L  of  water  was  decanted  from  the  sample  and  the  remaining 
2  L  of  activated  sludge  was  poured  into  Reactor  1.  A  manual  feed  was  initiated  using 
Feed  A  and  Feed  B  followed  by  the  start  of  an  automated  cycle. 

The  automated  cycle  covering  12  hrs  was  controlled  using  a  software  program, 
electronically  controlled  switches  and  solenoids  and  operated  continuously.  It  con¬ 
sisted  of  aeration,  settling,  effluent  purge,  influent  deionized  water  and  feed  solution 
addition.  During  the  10:50  hr  aeration  period,  filtered  compressed  air  was  pumped 
into  the  reactors  at  constant  pressure  to  maintain  approximately  7  ^  concentration 
dissolved  oxygen  (DO)  in  the  mixed  liquor.  Air  was  vented  out  during  the  aeration 
period,  except  during  the  effluent  purge,  which  occurred  5:30  hrs  into  the  cycle,  paus¬ 
ing  aeration  for  4  sec.  During  the  last  second,  the  35.9  mL  of  effluent  was  purged. 
At  the  end  of  10:50  hr  aeration  period,  the  mixed  liquor  was  allowed  to  settle  for 
1  hr.  After  settling,  a  purge  air  line  was  switched  on  and  the  supernatant  effluent 
was  pumped  for  45  sec  from  the  tanks  into  a  graduated  cylinder.  Then,  624  mL  of 
deionized  water,  36.0  mL  of  Feed  B,  and  8.0  mL  of  Feed  A  were  added  to  the  tank, 
finishing  at  the  end  of  the  12  hr  cycle.  Reactors  2  and  3  received  627  mL  of  deionized 
water  and  38.2  mL  of  Feed  B  during  operation  due  to  measured  delivery  differences 
noted  during  calibration.  Two  full  cycles  occurred  each  day  for  the  reactors  during 
the  operational  period. 

Reactor  1  activated  sludge  stabilization  and  adaptation  to  the  synthetic  influent 
continued  for  60  days.  The  wasted  mixed  liquor  purge  solution  was  collected  and 
added  to  Reactors  2  and  3  on  Day  70  of  operation.  Air  pressure  and  liquid  volumes 
were  calibrated  and  adjusted  to  account  for  unique  flow  rates  associated  with  each  re¬ 
actor  and  associated  lengths  of  tubing.  All  three  reactors  were  operated  continuously 
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at  room  temperature.  The  reactor  vessels  were  maintained  using  brushes  and  deion¬ 
ized  water  to  remove  biomass  material  from  the  sides  of  the  vessel.  Hoses  and  tubing 
were  periodically  inspected  and  cleaned  to  ensure  a  consistent  flow  rate.  A  Mettlcr 
Toledo  FiveEasy  bench  top  pH  meter  was  also  used  to  electrometrically  determine 
the  pH  of  the  mixed  liquor  solutions.  The  bench  top  meter  was  used  as  described  in 
Standard  Method  9040C,  pH  (Electrometric)  (26). 

After  Reactor  1  had  operated  for  250  days  and  all  three  reactors  had  stabilized 
using  the  synthetic  influent  feed  solution,  malathion  was  added  to  the  influent  feed 
for  Reactors  1  and  3.  Reactor  1  received  malathion-spiked  influent  at  a  rate  sufficient 
to  maintain  0.1^  concentration  and  Reactor  3  received  at  a  rate  to  maintain  3^A 
Reactor  2  did  not  receive  malathion  and  was  maintained  as  the  control.  The  influent 
water  volume  was  624  mL  per  cycle. 


mg  (0.1^)(2  L) 

Cone,  Malathion,  Rl,  — —  =  —  —  =  0.3205 

L  0.624L 

(0.3205^)  (4L)  1000mA 

Vol,  Malathion,  Rl,mL  =  - - nnmg - =  0.01282L( — — — )  =  12.82 


100^ 

mo  (3^)(2L) 

Cone ,  Malathion,  R3,  —  =  — - - 

’  ’  ’  L  0.624 L 


=  9.615 


(9.615^)(4L) 


IL 


,  lOOOmL 


Vol,  Malathion,  A3,  mL  =  - - nr^my0v — -  =  0.3846L( — — — )  =  384.6 
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3.3  COD  and  Nitrogen 


Chemical  Oxygen  Demand  (COD),  Nitrogen- Ammonia  (NH4-N),  Nitrogen-Nitrate 
(NO3-N),  and  Nitrite  (NO2-)  were  evaluated  using  Standard  Methods  8000,  10031, 
10020,  8153  (27),  respectively,  and  procedures  adopted  from  Janeczko  (1).  For  each 
method,  a  calibration  curve  consisting  of  five  measurements  at  known  concentrations 
was  generated  with  linear  regression  and  interpolation. 

Thirteen  each  COD  digestion  vials,  NH4-N  vials,  NO3-N  vials,  and  10  mL  test 
tubes  for  NO2-  tests  were  obtained  and  labeled  for  three  effluent  tests  and  one  mixed 
liquor  test  per  reactor,  plus  one  additional  tube  as  a  negative  control  for  all  three 
reactors.  Effluent  and  mixed  liquor  samples  were  obtained  from  the  gravimetric  anal¬ 
ysis  filtrate.  An  Agilent  Cary  60  UV-Vis  Spectrophotometer  was  used  to  measure 
absorbance  for  each  reaction.  The  absorbance  from  the  samples  was  identified  on  the 
calibration  curve  and  the  associated  concentration  recorded.  NO3-N  effluent  had  a 
dilution  of  2.  All  others  were  calculated  with  a  dilution  of  1.  A  single  mixed  liquor 
value  and  mean  effluent  value  for  the  three  COD,  NH4-N,  NO3-N,  and  NO2-  tests  for 
each  reactor  was  obtained. 

3.4  DNA  Extraction 

DNA  was  extracted  using  the  Mo  Bio  Laboratories,  Inc.  PowerSoil®  DNA  extrac¬ 
tion  kit  and  protocol  (28)  and  procedures  adopted  from  Janeczko  (1).  Nine  2  mL 
collection  tubes,  1  PowerBead  Tube,  and  1  Spin  Filter  Tube  were  obtained  for  each 
of  the  three  reactors,  labeled,  and  placed  in  a  microcentrifuge  tube  rack.  A  1  mL 
sample  was  obtained  from  the  unfiltered  25  mL  mixed  liquor  sample  and  centrifuged 
at  10,000  x  g  for  2  min  to  pellet  the  bacteria.  The  liquid  supernatant  was  then  re¬ 
moved  and  discarded.  The  buffer  solution  from  the  PowerBead  Tubes  was  then  added 
to  the  bacteria  pellet  and  vortexed  briefly  to  resuspend  the  pellet  in  solution.  The 
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mixed  buffer  and  pellet  were  then  added  back  to  the  PowerBead  tube.  60  jiL  of  Mo 
Bio  proprietary  disruption  agents  was  added  to  promote  cell  lysis  and  the  tube  was 
vortexed  briefly.  Each  tube  was  then  secured  horizontally  to  the  vortex  pad  with  tape 
and  vortexed  for  10  min  for  homogenization  and  cell  lysis. 

The  tube  was  then  centrifuged  at  10,000  x  g  for  30  sec  and  the  supernatant 
transferred  to  a  clean  collection  tube.  250  /iL  of  inhibitor  removal  reagents  was 
added  to  precipitate  contaminants  and  vortexed  for  5  sec  and  then  incubated  at  4°C 
for  5  min.  The  tube  was  then  centrifuged  at  10,000  x  g  for  1  min  and  the  supernatant 
transferred  to  a  clean  collection  tube.  200  /iL  of  additional  inhibitor  removal  reagents 
was  added  for  further  sample  purification  and  vortexed  for  5  sec  and  then  incubated 
at  4°C  for  5  min.  750  /xL  of  supernatant  was  transferred  to  a  clean  collection  tube 
and  1200  /i L  of  vortexed  high  concentration  salt  solution  was  added.  The  tube  was 
then  vortexed  for  5  sec.  675  /xL  of  the  solution  was  transferred  onto  a  spin  filter  with 
a  silica  membrane  and  centrifuged  at  10,000  x  g  for  1  min.  The  flow  through  was 
discarded  and  the  filter  transferred  to  a  clean  collection  tube. 

The  spin  filter  was  loaded,  centrifuged,  and  transferred  two  additional  times  with 
the  remaining  solution  for  a  total  of  3  sample  loads.  500  /xL  of  ethanol-based  wash 
solution  was  added  to  the  spin  filter  and  centrifuged  at  10,000  x  g  for  30  sec.  The  flow 
through  was  discarded  and  the  filter  transferred  to  a  clean  collection  tube.  The  spin 
filter  was  centrifuged  at  10,000  x  g  for  1  min  to  remove  all  residual  wash  solution  and 
transferred  to  a  clean  collection  tube.  100  /xL  of  sterile  elution  buffer  was  added  and 
centrifuged  at  10,000  x  g  for  30  sec.  The  spin  filter  was  removed  and  the  extracted 
sample  placed  in  -80° C  storage. 
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3.5  Sample  Quality 


The  extracted  samples  were  evaluated  for  contamination  using  a  Nanodrop™  1000. 
The  instrument  was  cleaned  with  deionized  water  prior  to  use.  A  blanking  buffer  of 
deionized  water  was  used  as  a  reference.  A  2  /j,L  drop  from  each  sample  was  pipetted 
onto  the  measurement  pedestal  and  the  sampling  arm  was  lowered  to  a  10  mm  path 
length.  A  measurement  was  taken  for  each  sample  and  the  pedestal  cleaned  with  a 
lab  tissue  afterwards.  Each  measurement  was  reviewed,  saved,  and  exported.  The 
A260/A280  and  A260/A230  ratios  were  examined  for  appropriate  values. 

Contaminants  from  reagents  used  in  sample  DNA  extraction  may  be  identified  by 
shifts  in  the  spectra  at  the  A260  peak,  A230  trough,  and  A280  value.  Examples  of 
contaminants  include:  TRIzol,  phenol,  guanidine  ffCL,  and  guanidinium  isocyanate. 
Low  A260/A230  ratios  may  indicate  residual  phenols,  guanidine,  or  glycogen,  ffigh 
A260/A230  ratios  may  indicate  inappropriate  blank  measurements.  Low  A260/A280 
ratios  may  indicate  residual  reagent  contamination,  such  as  phenols.  The  generally 
accepted  value  for  pure  DNA  is  1.8  for  A260/A280  and  slightly  higher  values  in  the 
range  of  2.0  -  2.2  for  A260/A230  (29). 

Three  samples  were  selected  from  Reactor  1  on  Days  105,  112,  and  119  of  operation 
(Fig.  2).  The  A260/A280  values  for  all  three  samples  ranged  between  1.83  and  2.07. 
The  A260/A230  values  ranged  between  0.57  and  1.7.  However,  the  low  A260/A230 
values  were  likely  affected  by  low  concentrations  of  dsDNA  in  the  samples  and  the 
proximity  to  the  lower  limits  of  detection.  The  detection  limits  for  the  Nanodrop™ 
1000  are  2  ng//xL  to  75  ng/yuL  at  the  10mm  pathlength  (30).  This  purity  values  at  the 
A260/A280  levels  and  the  location  of  the  peaks  were  within  the  proper  ranges.  High 
purity  DNA  was  obtained  using  the  previously  described  DNA  extraction  methods 
and  suitable  for  downstream  processing. 
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Figure  2.  Contaminant  Analysis 


3.6  Sample  Quantity 


Table  1.  PCR  product  concentrations  (before  malathion) 


Primer  Set 

Reactor  1  (ng/pL) 

Reactor  2  (ng/pL) 

Reactor  3  (ng/pL) 

Min,  Max,  Mean 

Min,  Max,  Mean 

Min,  Max,  Mean 

BaclLEXMl  /  BacU_EX_E786F 

8.070,  72.7,  43.55 

1.040,  68,  34.37 

0.887,  49.3,  27.38 

Bac-Ul_13  /  BaclLUX_U1053F 

31.50,  69.3,  50.91 

15.7,  63.1,  45.67 

14.1,  70,  39.17 

805r  /  Bact_U2_15 

1.10,  84,  26.65 

3.37,  175,  58.67 

4.94,  169,  65.55 

The  reactor  samples  were  evaluated  for  initial  quantity  using  the  Qubit®  2.0  flu- 
orometer  and  a  dsDNA  HS  (High  Sensitivity)  Assay  Kit.  Two  thin-wall,  clear  0.5 
mL  PCR  tubes  were  labelled  and  filled  with  189  /jL  of  Qubit™  working  solution  and 
1  /iL  of  Qubit™  dsDNA  HS  reagent.  10  /iL  of  0  ng//iL  Standard  #1  was  added  to 
one  tube  and  10  /iL  of  10  rig//iL  Standard  #2  was  added  to  the  other.  Both  were 
allowed  to  incubate  at  room  temperature  for  2  min  and  then  measured  to  calibrate 
the  instrument.  194  /j,L  of  working  solution,  1  /j,L  of  dsDNA  HS  reagent,  and  5  /i L  of 
each  sample  were  added  to  0.5  mL  PCR  tubes  and  incubated  at  room  temperature 
for  2  min.  Each  was  then  measured,  noted,  and  the  data  exported  to  a  hie. 

The  huorometer  uses  dyes  that  are  selective  and  only  bind  to  DNA,  minimizing 
the  effect  of  contaminants  on  the  quantitative  measurement  (31).  The  three  DNA 
samples  selected  for  analysis  before  PCR  ranged  between  0.133  and  16.36  ng///L 
(Fig.  3).  The  measurement  was  repeated  after  PCR  using  the  initial  primers  without 
barcodes  and  adapters  (Tab.  1).  For  the  medium  quantity  assay  (Day  113),  there 
was  a  29  -  77  times  increase  in  nucleic  material.  The  largest  quantity  assay  (Day  119) 
increased  by  9  -  30  times. 
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Qubit  Analysis 
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Figure  3.  Quantitation  of  Sample  Nucleic  Material 


3.7  Sequencing  Approach 


The  overall  objective  was  to  conduct  six  total  sequencing  iterations  representing 
samples  from  each  of  the  three  reactors  prior  to  malathion  exposure  and  after  expo¬ 
sure.  Broad  coverage  of  the  16S  ribosomal  gene  was  also  a  consideration.  A  multiplex 
approach  using  fragment  lengths  of  350  -  500bp  was  selected  initially.  Multiplex  se¬ 
quencing  was  attempted  on  the  first  sequencing  run,  but  later  runs  were  conducted 
using  a  single  primer  pair.  Early  evaluation  of  primer  performance  and  use  was  ac¬ 
complished  using  primers  with  manually  ligated  sequencing  adapters  and  barcodes. 
The  follow-on  evaluation  of  samples  from  before  and  after  malathion  exposure  was 
conducted  using  fusion  primers  with  sequencing  primers  and  barcodes  attached.  The 
barcoding  methodology  highlighted  evolution  of  the  bacterial  community.  Barcode  1 
represented  the  initial  period  of  operation  during  the  first  90  days.  Barcodes  2  and 
3  represented  the  early  and  late  periods,  respectively,  of  the  acclimation  period  prior 
to  malathion  exposure  on  Day  250.  Barcode  5  represented  the  30  days  of  malathion 
exposure.  Due  to  sampling  constraints,  Barcode  4  was  eliminated  from  the  final 
evaluation. 

3.8  Primer  Preparation  and  Use 

A  variety  of  primers  were  surveyed  for  use  during  sequencing.  Primers  were  se¬ 
lected  from  known  high  coverage  rate  pairs  from  previous  research  (32)  (33)  (34)  (35)  (36). 
An  example  bacteria,  Escherichia  coli,  NC_0 11 745.1,  was  examined  using  a  BLAST 
search  (37).  The  16S  ribosomal  rRNA  gene,  ECED1_16S_1  [3868436. .3870438],  was 
located  and  downloaded  using  the  graphical  interface.  The  sequence  was  used  to 
conduct  a  BLAST  search  versus  the  list  in  FASTA  format. 

The  offsets  and  strandedness  were  associated  and  mapped  on  a  linear  scale  repre¬ 
senting  the  16S  ribosomal  rRNA  gene  (Fig.  4).  Primer  pairs  were  reviewed  using  the 
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Table  2.  Primer  detail 


Sequence  (5’  to  3’) 

Sequence  Name 

Offset 

Strandedness 

RGYTACCTTGTTACGACTT 

BacU_EX_ll 

531 

+ 

GTGCCAGCMGCCGCGGTAA 

BacU_EX_E786F 

905 

- 

GAGTTTGATCMTGGCTCAG 

Bac-Ul_13 

217 

+ 

CCTACGGGNGGCWGCAG 

BacU_UX_U1053F 

640 

- 

ACTCCTACGGGAGGCAGCAGT 

805r 

903 

+ 

ATCGGYTACCTTGTTACGACTTC 

Bact_U2_15 

1368 

- 

Primer  Mix  2b 


52.7  deg  C 


217 

Bacll  EX 

531 

11 - 

594 

634 

640 

--374b  p- 
723 

767 

768 

902 

Bacll 

903 

X  E 786F 1 

905  1004 

1171 

1175 

1350 

1367 

1368 

1588 

1681 

1948 

♦ 

♦ 

♦ 

+ 

+ 

-  1 

♦ 

♦ 

Primer  Mix  3 


47.5  deg  C 


Bac-Ul  13 - 423bp ----- 

Bacll  UX 

U1053F 

217 

531 

594 

634 

640 

723 

767 

768 

902 

903 

905 

1004 

1171 

1175 

1350 

1367 

1368 

1588 

1681 

1948 

♦ 

+ 

+ 

+ 

+ 

♦ 

+ 

Primer  Mix  S 

_ 50.6  deg  C 


217 

531 

594 

634 

640 

723 

767 

768 

902 

805r - 

903 

905  1  1004 

1171  1 1175 

1350 

U2  15 (B 

1367 

act340F) 

1368 

1588 

1681 

1948 

♦ 

+ 

+ 

+ 

♦ 

'  1 

♦ 

+ 

Figure  4.  16S  rRNA  Primer  Alignment 


SILVA  rRNA  gene  database  for  bacterial  domain  coverage  (38).  Resulting  sequences 
were  then  filtered  using  fragment  lengths  of  350  -  500  bases  for  sequencing  on  the  Life 
Technologies  Ion  Torrent  Personal  Genome  Machine  (PGM)  using  the  316v2  chip. 
The  initial  BLAST  search  results  identified  alignments,  offsets,  and  strandedness  for 
the  list  of  primers  (Tab.  2). 

Initially,  six  printer  sets  were  selected  from  the  list  and  evaluated  using  E.  coli 
template.  Three  primer  pairs  with  overlapping  offsets  were  selected  based  upon  the 
quantity  of  PGR  product  and  length. 

3.9  PCR 

Non-barcoded,  lyophilized  primers  (200  nmol,  HPLC  purified)  were  obtained  from 
the  list  of  primer  pairs  identified  in  the  BLAST  alignment  (37)  and  16S  rRNA  ribo- 
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Table  3.  Gradient  PCR  Profile 


Step 

T  (°C) 

Time 

Cycles 

1 

50 

30min 

2 

94 

15min 

3 

94 

30sec 

4  (see  Fig.  5) 

47  -  65 

30sec 

5 

72 

lmin 

Repeat  3-5 

32x 

6 

72 

15min 

7 

4 

oo 

somal  gene  survey.  The  primers  were  hydrated  using  a  TE  Buffer  (10  niM  Tris»Cl/l 
mM  EDTA)  to  a  stock  solution  of  100  /iM  concentration.  Working  solutions  of  20 
/zM  were  prepared  from  the  stock  solution.  The  primer  solutions  were  then  placed  in 
-20°C  storage  until  use. 

Template  DNA  was  obtained  from  laboratory  stock  or  the  reactor  samples,  thawed, 
and  diluted  to  10  pg//iL  concentration.  PCR  reaction  mixtures  for  each  of  the  primer 
pairs  were  prepared  using  the  Qiagen  OneStep  RT-PCR  kit.  The  total  quantity  of 
reaction  mixture  was  calculated  for  an  eight  reactions.  120  /zL  DNA-grade  water,  50 
/zL  5X  buffer,  10  ytzL  10  mM  dNTP  mix,  5  ytzL  primer  A,  and  5  /zL  primer  B  were 
added  to  the  master  mix  tube.  Once  all  master  mixes  were  prepared,  10  /zL  of  en¬ 
zyme  was  added  to  each  tube  and  each  tube  was  kept  on  ice.  The  master  mixes  were 
then  aliquoted  into  eight  tube  strips  for  PCR.  Five  ytzL  of  template  was  added  to  each 
aliquot  for  a  total  of  25  ytzL  for  each  0.2  mL  tube  reaction. 


24 


CM 

PC? 

LO 

o 

00 


CO 

LO 

o 

CO 

T— I 
t-H 

2T 


<x> 

00 

i> 

w 


X 

H 


65.0 

X 

64.5 

X 

63.3 

X 

61.4 

X 

X 

59.3 

X 

X 

X 

57.1 

X 

X 

X 

54.9 

X 

X 

X 

52.7 

X 

X 

X 

50.6 

X 

X 

48.7 

X 

X 

47.5 

X 

X 

47.0 

X 

o 

H 

EX_11-E786F 

UL13-U1053F 

805R-U2.15 

CO 

CO 


fl 


fl 


T3 

cS 


u 


o 


u 

Ph 


in 

QJ 

Sh 

3 

be 


25 


A  temperature  gradient  was  used  to  identify  optimum  melting  temperatures  (Tm) 
for  the  initial  primers  without  adapters  or  barcodes  (Tab.  3).  Eight-strip  PCR  tubes 
with  prepared  template  DNA  and  reagent  were  placed  in  the  thermocycler.  The 
placement  in  the  gradient  PCR  wells  was  based  on  GC%  content  and  estimation 
from  previous  attempts  (Fig.  5). 

For  the  fusion  primers  with  adapters  and  barcodes,  the  PCR  temperature  profile 
was  determined  using  Integrated  DNA  Technologies’  OligoAnalyzer®  (39).  For  each 
search,  the  concentration  parameters  for  salts  and  dNTPs  were  updated  using  the 
Qiagen  OneStep  RT-PCR  kit  manual  as  a  reference.  The  primer  concentration  was 
input  as  0.25/zM.  The  primer  pair  sequences  only  (without  adapters  and  barcodes) 
were  input  into  the  query  box  and  the  resulting  Tm  were  reduced  by  7°C.  Then  the 
primer  pair  sequences  with  adapters  and  barcodes  were  input  and  resulting  Tm  were 
reduced  by  10°C.  These  values  were  used  in  PCR  for  the  fusion  primers  and  reactor 
sequencing  (Tab.  5).  The  PCR  profile  is  shown  in  Table  4. 

3.10  Gel  Electrophoresis 

The  products  from  PCR  were  evaluated  using  gel  electrophoresis.  The  tests  were 
run  on  E-Gel®  single  or  double  comb,  2%  agarose,  ethidium  bromide  stained,  precast 
gels  on  an  E-Gel®  PowerBase™  using  a  Traeklt™  50  bp  DNA  Ladder.  The  DNA 
ladder  was  diluted  1:4  in  a  1.5  mL  tube,  vortexed  and  placed  in  a  microcentrifuge 
tube  rack.  The  required  number  of  precast  gels  were  obtained  from  laboratory  stock 
and  indexed  on  the  bottom  right  corner  using  a  marker.  The  PCR  strips  containing 
template  DNA  were  removed  from  -80° C  storage  and  placed  on  ice.  For  the  single 
comb  tests,  lanes  2,  11,  and  12  were  filled  with  20  /zL  of  DNA  quality  water  and  lanes 
3-10  were  filled  with  10  /zL  of  DNA  quality  water.  Lanes  3-10  were  filled  with  10 
yuL  of  DNA  template.  Lane  1  was  filled  with  20  /iL  of  the  diluted  ladder. 
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Table  4.  Sequencing  Amplification  Profile 


Step 

T  (°C) 

Time 

Cycles 

1 

50 

30min 

2 

94 

15min 

3 

94 

30sec 

4 

T  1 

30sec 

5 

72 

lmin 

Repeat  3-5 

15x 

6 

94 

30sec 

7 

T  2 

J-  m 

30sec 

8 

72 

lmin 

Repeat  6-8 

25x 

9 

72 

15min 

10 

4 

oo 

Table  5.  Sequencing  Temperature  Determination 


Primer 

Tm1  (°C) 

Primer 

Tm2  (°C) 

Primer 

Adapter 

Barcode 

Bases  (max) 

GC  content  (%) 

BacU_EX_ll 

75 

81.6 

64 

54.7 

BacU_EX_E786F 

57.8 

78.3 

48 

54.2 

Selected  Tm(°C) 

50.8  (15x) 

68.3  (25x) 

Bac-Ul_13 

67.5 

79.4 

57 

51.8 

BacU_UXJJ1053F 

68.1 

81.6 

46 

63 

Selected  Tm(°C) 

60.5  (15x) 

69.4  (25x) 

805r 

63.1 

79.4 

56 

53.6 

Bact  _U2_15 

71.8 

81.7 

49 

61.2 

Selected  Tm(°C) 

56.1  (15x) 

69.4  (25x) 
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For  the  double  comb  tests,  lanes  2  -  4,  5  -  7,  10  -  12,  and  13  -  15  were  filled  with  10 
/jL  of  DNA  template  (Fig.  6).  Lanes  8  and  16  were  filled  with  E.  coli  template  as  the 
positive  control.  Lanes  1  and  9  were  the  negative  control.  Diluted  ladder  was  added 
to  the  middle  (M)  lane.  The  gel  was  then  placed  on  the  E-Gel®  PowerBase™  and  a 
30  min  timer  was  started.  After  30  min,  the  precast  gel  with  the  sample  was  removed 
and  placed  under  UV  transillumination  using  a  Kodak  Gel  Logic  200  Imaging  System. 
The  DNA  template  lanes  were  evaluated  for  clarity  at  the  target  lengths  and  optimal 
temperature  ranges  noted. 
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3.11  Adapter  and  Barcode  Ligation 


Once  successful  PCR  products  were  shown  using  gel  electrophoresis,  adapters  and 
barcodes  were  attached  to  the  DNA.  Ligation  was  only  used  when  fusion  primers  (with 
adapters  and  barcodes  already  attached)  were  not  used.  The  Ion  16S™  Metagenomics 
kit  workflow  was  used  as  guidance  for  workflow  for  the  remainder  of  the  library 
construction  (40).  The  Ion  Plus  gDNA  Fragment  Library  Preparation  kit,  barcode 
adapters  and  protocol  were  used  during  the  initial  phase  of  sequencing  parameter 
identification  (41).  The  reaction  was  setup  for  a  barcoded  library  using  additional 
polymerase.  41  fiL  of  nuclease-free  water  was  added  to  a  0.2  mL  PCR  tube.  10  /iL 
of  10X  ligase  buffer  was  added  to  the  reaction  mixture.  2  /iL  each  of  Ion  PI  adapter, 
Ion  Xpress™  Barcode  X,  and  dNTP  mix  were  then  added.  The  gDNA  input  was 
100  ng  at  25  yuL.  Finally,  2  yuL  of  DNA  ligase  and  16  yuL  of  nick  repair  polymerase 
were  added  bringing  the  total  reaction  volume  to  100  yuL.  The  reaction  tube  was  then 
place  in  a  thermocycler  programmed  with  a  preheated  lid  at  95°C  for  3  min,  25°C  for 
15  min,  72°C  for  5  min  with  a  final  hold  at  4°C.  The  reaction  was  removed  from  the 
thermocycler  once  the  temperature  reached  4°C  and  purified  using  the  Agencourt® 
AMPure®  XP  Reagent  kit. 

3.12  Purification 

Purification  using  the  Agencourt®  AMPure®  XP  Reagent  kit  consisted  of  attaching 
the  gDNA  to  beads,  removal  of  the  supernatant,  a  series  of  ethanol  washes,  and  elution 
using  a  TE  buffer.  Prior  to  beginning  the  purification  step,  a  fresh  mixture  of  70% 
ethanol  and  30%  nuclease-free  water  was  prepared.  The  sample  volume  was  measured 
using  a  pipette.  A  1.8X  sample  volume  (180  /xL)  of  AMPure®  beads  was  obtained 
from  4°C  storage  and  allowed  to  come  to  room  temperature  for  30  min  (41).  The 
sample  and  beads  were  combined  in  a  1.5  mL  Eppendorf  LoBind®  tube  and  mixed 
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using  a  200  n L  pipette  5  times.  The  mixture  was  incubated  at  room  temperature 
for  5  min.  The  mixture  was  pulse  spun  and  placed  in  a  magnetic  tube  rack  until  the 
solution  cleared.  The  supernatant  was  withdrawn  from  the  bead  pellet  and  500  rnL  of 
the  ethanol  mixture  was  added.  The  tube  was  rotated  slowly  in  the  rack  over  a  period 
of  30  sec  and  the  ethanol  wash  was  removed  using  a  pipette.  The  ethanol  wash  was 
repeated  a  second  time.  After  removal  of  the  wash  solution  using  a  pipette,  the  tube 
was  removed  from  the  rack,  pulse  spun,  a  placed  back  in  the  rack.  The  remaining 
ethanol  wash  was  removed  using  a  20  /jL  pipette.  The  tube  was  left  open  to  the 
air  for  a  period  of  4.5  min  to  allow  the  remaining  wash  to  evaporate.  The  tube  was 
removed  from  the  rack  and  40  /iL  of  TE  buffer  was  added  directly  to  the  pellet.  The 
buffer  and  beads  were  mixed  using  a  pipette  10  times.  The  tube  was  pulse  spun  and 
placed  in  the  magnetic  rack  until  the  solution  cleared.  The  supernatant  contained 
the  gDNA  was  removed  in  two  20  /iL  volumes  into  2  LoBind®  tubes.  Both  tubes  were 
placed  in  -80° C  storage. 

3.13  Size  Selection 

Size  selection  of  the  purified  and  adapter  ligated  and  barcoded  gDNA  was  con¬ 
ducted  using  the  Pippen  Prep®  instrument  and  2%  agarose  and  ethidium  bromide  gel 
cassettes  with  external  marker  B  and  manufacturer’s  protocol  (42).  Prior  to  opera¬ 
tion,  the  reagents  for  the  kit  were  removed  from  4°C  storage  and  allowed  to  come  to 
room  temperature  for  30  min.  Loading  buffer  was  added  to  the  sample  to  bring  the 
sample  volume  to  40  yuL.  The  Pippen  Prep®  instrument  was  calibrated  using  the  cal¬ 
ibration  optical  fixture.  The  cassette  was  removed  from  its  packaging  and  inspected 
for  proper  condition.  It  was  tilted  and  all  air  bubbles  were  dislodged  to  the  right  side 
of  the  cassette.  The  cassette  was  placed  in  the  instrument  tray  and  both  tape  tabs 
were  removed.  Approximately  40  -  50  /i L  was  removed  from  each  elution  well  and 
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replaced  with  fresh  electrophoresis  buffer.  The  tape  tab  from  the  kit  was  then  placed 
over  the  elution  wells.  Fresh  electrophoresis  buffer  was  added  to  each  sample  port  un¬ 
til  full.  A  current  test  was  performed  and  all  wells  of  the  cassette  passed.  A  program 
with  cut  sizes  matching  the  base  pair  size  targets  ±  35  bp  was  added  and  loaded  for 
the  cassette.  The  reference  lane  was  designated  as  lane  1  and  applied  to  all  sample 
lanes.  The  instrument  was  opened  and  40  /iL  removed  from  each  of  the  sample  lanes 
and  the  reference  lane.  40  / zL  of  sample  was  added  to  each  sample  lane  and  40  /zL 
of  reference  ladder  was  added  to  the  reference  lane.  The  instrument  lid  was  closed 
and  the  program  was  started  and  ran  for  approximately  1.5  hours.  Once  the  program 
completed,  the  sample  fractions  (about  40yuL)  were  removed  from  the  elution  wells. 
The  fractions  were  purified  using  the  Agencourt®  AMPure®  XP  Reagent  kit.  The 
purification  protocol  was  identical  to  the  previous  protocol  (41).  After  purification, 
the  samples  were  placed  in  -80°C  storage. 

3.14  Quality  Control  and  Normalization 

Quality  control  of  the  samples  was  conducted  using  the  Agilent  Technologies  Bio¬ 
analyzer  2100  and  protocol  (43).  Samples  were  removed  from  -80°C  storage  and, 
along  with  the  high  sensitivity  DNA  analysis  reagents,  were  allowed  to  come  to  room 
temperature  for  30  min  prior  to  beginning  the  analysis.  The  electrode  assembly  was 
removed  from  the  Bioanalyzer  and  sonicated  for  10  min  and  dried  using  compressed 
air.  The  electrode  assembly  was  returned  to  the  Bioanalyzer  system.  A  chip  was 
loaded  into  the  chip  priming  station  and  the  plunger  on  the  priming  syringe  was  set 
to  1  mL.  9  /iL  of  gel-dye  mix  from  the  high  sensitivity  DNA  analysis  reagent  kit 
was  added  to  the  gel  priming  well.  The  priming  station  was  closed  and  the  syringe 
depressed  until  locked  into  place  by  the  priming  station  clip.  Exactly  1  min  later,  the 
clip  was  released  and  5  sec  later  the  syringe  plunger  was  slowly  withdrawn  to  the  1 
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Table  6.  Region  Table 


Sample  Name 

Average  Size  [bp] 

Cone,  [pg/1] 

Molarity  [prnol/1] 

%  of  Total 

EX_11-E786F 

466 

30,269.88 

98,548.50 

40 

U1_13-U1053F 

395 

22,200.29 

85,826.00 

71 

805R-U2_15 

512 

1,226.50 

3,636.30 

53 

mL  position.  The  chip  priming  station  was  opened.  9  pL  of  gel-dye  mix  was  added  to 
the  remaining  gel  wells.  6  pL  of  marker  was  added  to  all  wells  without  samples  and 
5  pL  added  to  all  sample  wells.  1  /i L  of  ladder  reagent  was  added  to  the  ladder  well. 
Finally,  1  pL  of  each  sample  was  added  to  each  sample  well.  The  chip  was  placed 
in  a  chip-specific  vortex  at  2,000  rpm  for  1  min.  The  chip  was  immediately  loaded 
into  the  Bioanalyzer  2100  system.  The  Bioanalyzer  program  was  programmed  to  use 
the  high  sensitivity  DNA  analysis  protocol,  the  program  was  initiated  and  ran  for  45 
min. 

Using  the  analysis  from  the  quality  control  using  the  Bioanalyzer  2100,  the  gDNA 
input  from  each  of  the  the  primer  pair  PCR  products  was  normalized.  Each  of 
the  products,  after  size  selection  and  purification,  were  evaluated  according  to  the 
molarity  measurement  from  the  region  table  from  the  Bioanalyzer  instrument  (Tab. 
6).  The  region  analysis  was  conducted  by  assigning  ranges  based  upon  the  peak 
readings  at  the  highest  values  displayed  (Fig.  7).  The  low  value  before  the  rise  in  the 
highest  peak,  The  high  value  was  after  highest  peak  declined  to  a  value  close  to  the 
null  value.  The  region  quantitation  was  input  into  a  spreadsheet.  A  dilution  value 
was  assigned  to  normalize  each  of  the  values.  The  molarities  were  then  equalized 
to  a  volume  ±  1  /i L  at  10  pL  of  common  volume.  Each  primer  was  normalized  to 
equivalency  at  45  -  52  pM  concentration  for  library  input  into  emulsion  PCR. 
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3.15  Emulsion  PCR 


The  Ion  One  Touch™  instrument  was  used  to  amplify  the  solution  for  sequencing 
using  the  manufacturer’s  protocol  (44).  The  lid  of  the  instrument  was  opened  and 
a  recovery  tube  installed  into  the  slots.  A  recovery  router  was  installed.  The  instru¬ 
ment  lid  was  closed.  A  new  amplification  plate,  tubing,  and  disposable  injector  were 
installed.  Recovery  solution  and  oil  were  added  to  bring  the  volumes  to  the  desig¬ 
nated  levels.  All  reagents  were  allowed  to  come  to  room  temperature  prior  to  use. 
The  reagent  mix,  Reagent  B,  Enzyme  Mix,  and  Reagent  X  were  prepared  and  added 
to  the  solution  according  to  the  protocol.  25  /jL  of  normalized  solution  was  added 
to  the  mixture  for  a  total  volume  of  900  /i L.  The  mixture  was  combined  by  mixing 
using  a  pipette.  100  /iL  of  Ion  Sphere  Particles  (ISPs)  was  added  to  the  mixture.  The 
mixture  was  combined  by  mixing  using  a  pipette.  The  complete  solution  was  then 
added  to  a  Ion  PGM™  OneTouch™  Plus  Reaction  Filter  Assembly  via  the  sample 
port.  1500  /i L  of  reaction  oil  was  then  added  via  the  sample  port.  The  filter  assembly 
was  then  inverted  and  installed  on  the  OneTouch™  instrument.  The  OneTouch™ 
instrument  was  then  initiated  using  the  Ion  PGM™  Template  OT2  400  Kit  protocol 
assist  and  ran  for  7.5  hrs.  A  final  spin  was  conducted  at  the  completion  of  the  run 
and  the  ISPs  processed.  Processing  included  removal  of  the  recovery  tubes  and  solu¬ 
tion  to  100  /xL.  The  ISPs  from  both  recovery  tubes  were  suspended  in  500  /iL  of  Ion 
OneTouch™  Wash  Solution,  dispersed,  and  stored  at  4°C. 

3.16  ISP  Quality  Control  and  Enrichment 

The  ISPs  were  incubated  at  50° C  for  2  min.  The  ISPs  were  then  centrifuged  at 
15,500  x  g  for  2.5  min  and  all  but  100  yuL  of  wash  solution  removed.  The  quality  of 
the  the  ISPs  was  evaluated  using  the  Qubit®.  The  Ion  Sphere™  Quality  Control  kit 
Ion  Probes,  Annealing  Buffer,  and  Quality  Control  Wash  Buffer  were  thawed.  A  2 
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/xL  aliquot  of  the  unenriched  ISPs  was  tranferred  to  a  0.2  mL  Eppendorf  PCR  tube 
and  brought  to  room  temperature  with  98  /iL  of  Ion  OneTouch™  Wash  Solution.  2 
/xL  was  transferred  to  a  0.2  mL  Eppendorf  PCR  tube.  19  /xL  of  annealing  buffer,  1 
/iL  of  Ion  Probes  were  added  to  the  tube  and  mixed  with  a  pipette.  The  tube  was 
loaded  into  a  thermocycler  with  a  heated  lid  set  a  95°C,  95°C  for  2  min,  and  37°C  for 
2  min.  The  unbound  probes  were  washed  three  times  with  200  /xL  of  Quality  Control 
Buffer  by  adding  200  /xL,  mixing  with  a  pipette,  centrifuging  at  15,500  for  1.5  min  and 
removing  the  supernatant.  10  /xL  was  left  in  the  tube  using  a  comparison  tube  as  a 
reference.  After  the  final  wash,  the  volume  was  brought  to  200  /xL  and  mixed  using  a 
pipette.  The  entire  volume  was  tranferred  to  a  Qubit®  assay  tube.  A  negative  control 
was  produced  using  a  similar  volume  without  sample.  Alexa  Fluor®  standards  488 
and  647  were  measured  on  the  Qubit®  fluorometer.  The  negative  control  and  sample 
containing  the  ISPs  were  then  measured  in  both  488  and  647  settings.  The  values 
for  the  standards,  control,  and  sample  were  then  input  into  the  calibration  factor 
spreadsheet  from  the  manufacturer’s  website  to  calculate  the  percent  templated  ISPs. 
A  value  between  10  -  30%  was  the  optimal  templated  ISPs  goal.  The  manufacturer’s 
recommended  molarity  of  26  pM  appeared  to  correlate  with  the  lower  end  of  this 
optimal  range.  A  molarity  of  26  -  52  pM  was  utilized  during  this  study. 

Following  a  successful  quality  control  evaluation,  the  ISPs  were  enriched  on  the 
Ion  OneTouch™  Enrichment  System  (ES)  using  the  manufacturer’s  protocol  (44).  A 
1  mL  volume  of  1M  NaOH  was  prepared  by  diluting  10M  NaOH  using  water  from  the 
sequencing  kit.  A  second  1  mL  volume  of  lOOmM  NaOH  was  prepared  in  the  same 
way  and  set  aside  for  use  later  in  sequencing.  A  melt-off  solution  was  prepared  by 
combining  280  /xL  of  Tween®  solution  with  40  /jL  of  1M  NaOH.  Dynabeads®  MyOne™ 
Streptavidin  Cl  Beads  were  removed  from  storage  and  resuspended  using  a  vortex. 
13  yuL  were  removed  and  added  to  a  1.5  mL  Eppendorf  LoBind®  tube.  The  tube 
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was  then  placed  on  a  magnetic  rack  for  2  min  and  the  supernatant  removed.  130  /zL 
of  MyOne™  Bead  Wash  Solution  was  added,  vortexed  and  briefly  centrifuged.  The 
8-well  strip  from  the  Ion  OneTouch™  kit  was  prepared  with  the  square-shaped  tab  on 
the  left.  Wells  3,  4,  and  5  were  filled  with  300  /j>L  of  Ion  OneTouch™  wash  solution. 
Well  7  was  filled  with  320  /iL  of  the  freshly  prepared  melt-off  solution.  Well  2  was 
filled  with  130  /jL  of  MyOne™  Beads  Wash  Solution  with  MyOne™  Beads.  Well  1 
was  filled  with  the  entire  volume  of  sample  template-positive  ISPs  (~100  /iL)  from 
the  emulsion  PCR  step.  The  strip  was  loaded  into  the  right  side  of  the  slot  in  the  tray 
of  the  enrichment  system  with  the  square  tab  to  the  left.  A  new  tip  was  installed  in 
the  tip  arm  and  a  0.2  mL  with  10  ytzL  of  neutralization  solution  placed  in  the  hole  of 
the  tip  loader.  The  enrichment  cycle  was  initiated  and  ran  for  approximately  35  min. 
Immediately  following  enrichment,  the  solution  was  mixed  in  the  PCR  tube  with  the 
neutralization  solution  by  pipetting  up  and  down  10  times.  The  PCR  tube  with  the 
enriched  beads  was  then  placed  on  ice. 

3.17  Sequencing 

The  sequencing  solutions  for  Wash  2  and  Wash  3  and  annealing  buffer  were  re¬ 
moved  from  storage  and  brought  to  room  temperature  in  a  closed  drawer.  The  dNTPs, 
sequencing  primer,  and  control  Ion  Sphere™  Particles  (ISPs)  were  placed  in  an  ice 
bucket  to  thaw.  A  cleaning  cycle  was  performed  on  the  Ion  PGM™  sequencer  using 
either  the  chlorite  protocol  or  the  the  water  rinse.  The  chlorite  protocol  was  used  once 
per  week  and  the  rinse  used  for  daily  uses.  A  previously  designated,  used  initialization 
Ion  316™  v2  chip  was  placed  in  the  sequencer  for  the  initialization  step.  The  three 
wash  bottles  from  the  400  kit  were  washed  3  times  using  18.212  water  and  inverted  to 
dry.  The  wash  bottle  sipper  tubes  were  then  removed  and  the  400  kit  selected  on  the 
sequencer  screen.  The  sequencer  verified  gas  pressure  after  screen  confirmation  and 
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a  glove  change  was  completed.  New  sipper  tubes  with  the  gray  tops  were  attached 
with  the  long  W2  tube  resting  on  a  sterile  gauze  patch.  350  /iL  of  100  mM  NaOH 
was  added  to  the  Wash  1  bottle,  which  was  then  capped.  50  mL  of  Wash  3  solution 
from  the  kit  was  added  to  the  Wash  3  bottle  and  the  bottle  capped.  2  L  of  18.211 
water,  the  entire  bottle  of  W2  solution  from  the  kit,  and  70  / iL  of  100  mM  NaOH  was 
added  to  the  Wash  2  bottle,  which  was  then  capped  and  mixed  by  inverting  5  times. 
The  Wash  2  bottle  was  then  attached,  followed  by  the  Wash  1  and  Wash  3  bottles 
and  the  caps  secured.  The  initialization  process  continued.  The  pH  was  monitored 
until  initialization  was  complete  to  allow  for  manual  adjustment  using  titration. 

After  initialization  completed  successfully,  the  next  step  was  preparing  the  dNTP 
bottles.  The  dNTP  reagent  bottles  were  designated  using  their  respective  stickers. 
20  /iL  of  each  dNTP  was  added  to  its  designated  bottle  by  placing  the  pipette  tip 
approximately  one  third  of  the  way  down  against  the  inside  wall  and  depressing 
the  pipette  plunger.  The  bottles  were  then  capped  and  placed  in  racks.  The  old 
sipper  tubes  were  removed  from  the  dNTP  ports  and  a  glove  change  was  conducted. 
New  sipper  tubes  with  blue  tops  and  the  appropriate  dNTP  bottles  were  then  firmly 
secured  to  the  ports.  The  machine  then  conducted  a  leak  test  and  the  dNTP  bottles 
filled. 

After  the  dNTP  bottles  were  prepared,  the  control  Ion  Sphere™  particles  were 
added  and  the  sequencing  primer  was  annealed.  The  control  particles  were  vortexed 
and  centrifuged  and  5  /iL  of  control  particles  added  to  the  PCR  tube  from  the  enrich¬ 
ment.  The  solution  was  then  mixed  with  a  pipette  and  centrifuged  for  2  min  at  15,500 
x  g.  The  supernatant  was  removed,  leaving  15  /iL  in  the  tube  by  using  another  tube 
for  comparison.  12  /i L  of  sequencing  primer  was  added.  The  total  volume  was  verified 
using  a  pipette  and  additional  annealing  buffer  added  as  required.  The  solution  was 
then  mixed  using  a  pipette.  The  PCR  tube  was  then  loaded  into  a  thermocycler  at 
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95° C  for  2  min  and  37° C  for  2  min  using  the  heated  lid  option. 

While  the  PCR  program  was  running,  the  sequencing  run  was  programmed  using 
the  administrative  interface  in  a  web  browser.  The  appropriate  chip  type  (316v2), 
generic  sequence,  and  no  reference  library  were  selected.  The  run  and  sample  iden¬ 
tifications  were  input,  along  with  designating  the  appropriate  kits  used.  Finally,  the 
program  was  saved  and  verified  under  the  “Plan  run”  selection. 

The  last  preparation  step  prior  to  initiating  the  sequencing  run  was  to  conduct  a 
chip  check  and  load  the  chip.  At  this  point,  gloves  were  removed  to  prevent  static 
discharge  that  would  damage  the  chip.  The  initialization  chip  was  removed  from  the 
sequencer  after  touching  the  grounding  plate.  A  new  316v2  chip  was  removed  from 
its  packaging  and  placed  in  the  sequencer.  The  barcode  was  read  by  the  scanner 
and  chip  check  initiated  on  the  touchscreen.  After  a  successful  chip  check,  sequencing 
polymerase  was  removed  from  storage.  The  sample  with  ISPs  and  primer  was  removed 
from  the  thermocycler.  3  /iL  of  polymerase  was  added  to  the  sample  by  rotating  the 
pipette  pusher  to  avoid  any  bubbles.  The  sample  with  polymerase  was  then  mixed 
using  the  pipette  and  allowed  to  incubate  at  room  temperature  for  5  min.  The 
sequencing  chip  was  then  removed  from  the  sequencer  and  all  liquid  removed  from 
the  chip.  The  loading  procedure  used  was  the  alternate  method  using  simplified 
instructions  and  the  weighted  chip  bucket.  This  involved  centrifuging  the  chip  upside 
down  in  the  weighted  bucket  for  5  sec.  The  30  /i L  sample  was  then  loaded  into  the 
chip  by  rotating  the  pipette  pusher  at  a  rate  of  1  /iL  per  second  to  prevent  any  bubbles 
entering  the  chip.  The  chip  was  then  placed  in  the  weighted  bucket  and  centrifuged 
with  the  chip  tab  in  for  1  min  and  then  centrifuged  again  with  the  chip  tab  out  for  1 
min.  The  chip  tab  was  then  tapped  on  the  bench  top  3  times  and  liquid  was  removed 
from  the  chip  loading  port  by  rotating  a  pipette  pusher  until  no  more  liquid  entered 
the  pipette  tip. 
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The  sequencing  run  was  then  initiated  by  selecting  the  run,  loading  the  chip, 
and  initiating  the  run.  The  planned  run  that  had  been  previously  programmed  was 
selected,  the  settings  reviewed,  and  confirmed.  After  grounding,  the  chip  was  loaded, 
clamped  and  the  cover  closed.  The  sequencer  calibrated  the  chip  and  the  sequencing 
run  was  initiated.  The  sequencing  run  time  was  approximately  5  hours.  The  process 
of  constructing  a  library,  preparing  the  templated  ISPs  and  sequencing  was  repeated 
for  each  reactor  two  times. 

3.18  Bioinformatics 

After  sequencing,  a  Binary  Alignment/Map  (BAM)  file  was  generated  with  raw 
sequencing  data  from  the  run.  The  SAMtools  0.1.18  utilities  suite  was  used  to  interact 
with  and  conduct  immediate  post-run  processing  of  the  BAM  file  (45).  Sequence 
adapter  trimming  was  accomplished  using  Cutadapt  1.2.1  (46).  Quality  control  and 
assembly  of  the  sequence  reads  was  conducted  using  Mira  3. 9. 1.7  (47).  The  resulting 
three  output  files  were  generated  with  FASTA,  FASTQ,  and  QUAL  extensions. 

A  standalone  version  of  NCBI’s  BLAST  2.2.30+  was  downloaded  along  with  the 
16S  Microbial  database  (37).  The  add-on  taxonomy  database  was  also  loaded  to 
include  additional  information  for  reference.  Using  the  command  line  interface,  the 
FASTA  files  for  each  barcode  set  of  contiguous  (contig)  sequences  were  parsed  using 
a  nucleotide  query  (blastn)  and  NCBI’s  16S  microbial  database.  The  BLAST  query 
used  default  standalone  values  and  was  output  in  XML  format. 

The  BLAST  output  file  was  then  imported  into  MEGAN  (48).  A  Least  Common 
Ancestor  (LCA)  parameter  filter  was  applied.  The  minimum  bit  score  for  results  was 
50.0.  The  maximum  expected  value  (e-value)  was  0.01.  All  returned  values  below 
10.0%  of  the  maximum  percentage  hit  were  discarded.  No  minimum  percentage  or 
number  of  reads  was  set.  Finally,  the  16S  Percent  Identity  Filter  was  enabled  enforcing 
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rank-based  match  percentage  requirements  of  Species  99%,  Genus  97%,  Order  90%, 
Class  85%,  and  Phylum  80%.  Under  the  Taxonomy  tab,  the  Parse  Taxon  Names 
feature  was  enabled.  The  BLAST  output  hie  was  then  imported.  Once  all  BLAST 
outputs  were  imported  for  a  reactor,  the  resulting  trees  were  compared  to  assess  the 
differences  between  reactor  taxa  compositions  and  supporting  contig  counts. 
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IV.  Results  and  Analysis 


4.1  Sequencing 


Table  7.  Sequencing  Run  Data 


Reactor  1 

Reactor  2 

Reactor  3 

ISP  Loading 

59% 

72% 

56% 

Total  Bases  Read 

116,000,000 

560,000,000 

404,000,000 

ISP  Not  Templated 

(10%) 

0% 

(8%) 

ISP  Polyclonal 

(36%) 

(25%) 

(36%) 

ISP  Low  Quality 

(66%) 

(8%) 

(20%) 

Useable  Reads 

728,779 

3,128,271 

1,680,002 

Useable  Reads  (%) 

22% 

69% 

52% 

Read  Mean  Length/Median/Mode  (bp) 

160  /  88  /  29 

179  /  99  /  41 

241  /  242  /  450 

Total  Contigs 

497 

484 

233 

Mean  Contig  Length  (bp) 

416 

395 

385 

Mean  Reads  per  Contig 

338 

326 

392 

Six  sequencing  runs  were  conducted  in  total.  Three  runs  were  selected  after  fil¬ 
tering  by  post-sequencing  quality  checks  (Tab.  7).  The  initial  total  number  of  reads 
was  filtered  to  ensure  the  optimal  quality  for  downstream  analysis.  The  most  efficient 
product  was  obtained  from  Reactor  2,  producing  69%  useable  reads  from  72%  ISP 
loading  and  560  million  reads.  Reactor  1  produced  the  most  contiguous  sequences 
at  497  total  contigs  after  sequencer-  and  software-based  quality  checks.  Reactor  3 
produced  the  fewest  contigs  at  233,  but  had  the  most  overall  mean  reads  per  contig 
at  392,  and  the  highest  read  length  parameters.  A  graphical  summary  is  shown  in 
Figures  8(a)  -  10(a).  The  highest  ISP  density  is  shown  in  red,  with  decreasing  values 
in  yellow,  green  and  blue,  respectively.  The  proportion  of  the  run  length  associated 
with  the  target  ~490bp  region  is  identified  by  a  peak  on  the  right  of  the  histogram 
(Fig.  8(c)  -  10(c)). 
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(a)  ISP  Density 

Figure  8. 


(a)  ISP  Density 

Figure  9. 


(a)  ISP  Density 

Figure  10. 


(b)  ISP  Summary 


(c)  Read  Length 


Reactor  1  Sequencing  Data 


(b)  ISP  Summary 


(c)  Read  Length 


Reactor  2  Sequencing  Data 


56% 

Loading 


3,553,609 


44% 

ipty  Wells 


92% 

Enrichment 


64% 


3,253,844 

2,091,041 


8% 

|No  Template 


36% 

Polyclonal 


0%  Test  Fragments 
0%  Adapter  Dimer 
20%  Low  Quality 


(b)  ISP  Summary 


(c)  Read  Length 


Reactor  3  Sequencing  Data 
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4.2  Nitrification 


Reactor  nitrification  and  suspended  solids  data  are  summarized  in  Figures  11  -  13. 
Average  values  for  each  were  compared  for  the  early,  acclimation,  late,  and  malathion 
addition  periods  between  Days  1  -  90,  91  -  170,  171  -  250,  and  251  -  283,  respectively. 
The  average  sludge  retention  time  (SRT)  for  across  all  periods  for  all  three  reactors 
was  19  (±6)  days.  The  average  pH  values  for  Reactors  1,  2,  and  3  were  7.16,  6.99, 
and  7.27,  respectively. 

In  all  three  reactors,  the  ratio  of  effluent  TSS  to  mixed  liquor  TSS  was  high  during 
the  early  period.  In  Reactors  1  and  3,  this  coincided  with  high  values  of  TSS  in  the 
effluent  at  124  ^  and  70  respectively,  and  in  Reactor  2  with  low  values  for  TSS 
in  the  mixed  liquor  at  1144  NCH-  concentrations  were  high  for  all  reactors,  but 
more  so  in  Reactors  2  and  3,  with  82%  of  nitrite  remaining  in  the  effluent  for  both 
compared  to  the  starting  mixed  liquor  value  at  the  beginning  of  a  12-hr  cycle.  This 
may  have  indicated  low  overall  proportions  of  NOB,  with  the  effects  amplified  by 
subsampling  to  seed  Reactors  2  and  3.  Finally,  COD  concentrations  in  the  effluent 
were  highest  for  Reactor  3  for  this  period  at  23.99  suggesting  a  lower  proportion 
of  heterotrophs.  The  results  show  a  disruption  of  the  recently  established  equilibrium 
achieved  in  Reactor  1. 

During  the  acclimation  period,  the  total  biomass  in  Reactor  1  dropped  to  74%  of 
the  higher  mixed  liquor  TSS  value.  In  Reactor  2,  a  higher  than  average  value  of  50% 
of  COD  remaining  the  effluent  compared  to  the  mixed  liquor  was  observed.  NH4-N 
concentrations  were  low  in  the  reactor  and  zero  in  the  effluent.  Nitrate  levels  were 
lower  in  the  effluent  than  in  the  mixed  liquor,  which  may  have  been  the  result  of 
denitrification  activity.  In  Reactor  3,  lower  than  average  NH4-N  removal  and  excess 
NO2-  in  the  mixed  liquor  and  effluent  indicated  degraded  nitrifying  activity.  These 
measurements  show  that  the  seeding  of  Reactor  2  and  3  involved  unequal  distribution 
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of  nitrifying  and  heterotrophic  groups  into  those  reactors. 

In  the  late  period  just  before  the  addition  of  malathion,  SRT  was  lower  in  all  three 
reactors  at  14  days.  NH4-N  and  NO2-  levels  in  Reactor  1  reached  a  peak  value  in 
the  starting  mixed  liquor  value  in  comparison  to  the  other  periods,  but  the  values  of 
both  were  zero  in  the  effluent.  This  corresponded  with  an  average  of  29  ^  nitrate 
in  the  effluent.  In  Reactor  2,  all  values  were  within  one  deviation  of  the  mean  across 
all  periods.  In  Reactor  3,  the  mean  effluent  COD  concentration  dropped  to  its  lowest 
value  at  15.86  indicating  an  active  heterotroph  community.  Effluent  NH4-N  and 
NO2-  concentrations  were  both  zero  and  nitrate  levels  averaged  28.76  For  all 
three  reactors,  the  communities  appeared  to  have  reached  an  equilibrium  structure. 

Finally,  malathion  was  added  to  the  feed  of  Reactor  3  on  Day  242  with  a  con¬ 
centration  that  achieved  3.0  in  the  mixed  liquor.  In  Reactor  1,  it  was  added  to 
achieve  0.1  ^  on  Day  255.  Monitoring  continued  until  Day  283.  In  Reactor  1,  NO2- 
concentrations  were  present  in  the  effluent  6  days  after  malathion  was  added,  but 
returned  to  zero  5  days  later.  Average  effluent  COD  levels  for  this  time  period  were 
lower  at  12.5  These  measurements  indicate  some  level  of  negative  impact  to  the 
NOB,  but  continued  metabolic  activity  among  the  heterotrophs.  However,  Reactor  2, 
which  was  maintained  without  adding  malathion,  showed  low  effluent  COD  concen¬ 
trations  and  expected  concentrations  of  0  ^  for  both  NH4-N  and  NO2-.  In  Reactor 
3,  the  mixed  liquor  COD  concentration  was  elevated  at  65  It  may  be  that  at 
3.0  ^  malathion,  heterotrophs  begin  to  exhibit  a  stress  response  towards  the  new 
substrate  as  described  in  Janeczko  (2014)  (1). 
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Effluent  DAY  255 
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Figure  11.  Reactor  1  Nitrogen,  Suspended  Solids,  and  pH 


Fffluent 
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Figure  12.  Reactor  2  Nitrogen,  Suspended  Solids,  and  pH 


Effluent  DAY  242 
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Figure  13.  Reactor  3  Nitrogen,  Suspended  Solids,  and  pH 


4.3  Bioinformatics 


Table  8.  Phylogenetic  Analysis  Results 


Reactor  1 

Reactor  2 

Reactor  3 

Early/Acclim/Late  or  Malathion  (Diversity) 

15.3  /  11.6  /  8.0 

9.1  /  4.5  /  3.2 

5.0  /  4.0  /  3.5 

Genera 

37 

27 

22 

Genera  (Reads  /  %  of  total) 

67  (32%) 

65  (36%) 

50  (40%) 

The  bioinformatics  analysis  of  the  three  reactors  described  the  similarities  and 
differences  between  the  reactors  and  the  early,  acclimation,  and  late  or  malathion 
addition  periods.  The  biodiversity  of  each  sample  set  was  assessed  using  the  Simpson 
Reciprocal  Index  feature  in  MEGAN  (48).  The  most  diverse  sampling  was  obtained 
from  Reactor  1  with  values  of  15.364,  11.636,  and  8.000  for  the  acclimation,  malathion 
addition,  and  early  periods,  respectively  (Tab.  8).  The  least  diverse  sampling  was 
from  Reactor  3  with  values  of  5.000,  4.000,  and  3.571,  from  the  acclimation,  malathion 
addition,  and  early  periods,  respectively.  Reactor  2  sample  values  were  9.151,  4.500, 
and  3.267  from  the  acclimation,  early,  and  late  periods.  These  values  indicate  the 
probability  of  obtaining  duplicate  genera  in  any  random  sampling  from  the  reactors 
ranged  between  6.5%  to  30.61%.  All  three  reactors  were  associated  with  a  temporal 
decline  in  diversity  during  the  operational  period. 

For  each  period  and  reactor,  a  metagenomic  comparison  was  conducted  in  MEGAN 
to  associate  reads  with  taxa.  MEGAN  uses  similarity-based  binning  and  reference  se¬ 
quences  from  known  operational  taxonomic  units  (OTUs)  to  build  a  tree.  Sequences 
without  annotation  are  placed  between  annotated  sequences,  and  based  on  similarity, 
the  unannotated  sequence  is  placed  nearer  to  the  more  similar  annotated  sequence. 
Assignments  of  reads  to  taxons  is  then  conducted  using  LCA  deduction,  with  non¬ 
homologs  being  assigned  to  the  “no  hits”  category  (49). 

The  trees  shown  in  Figures  14  -  16  display  taxonomic  profiles  across  all  three 
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periods.  A  bar  chart  next  to  each  taxon  displays  the  supporting  square  root  nor¬ 
malization  of  reads,  a  recommended  statistical  method  for  phylogenetic  analysis  (50). 
The  bars  represent  the  early,  acclimation,  and  late/malathion  addition  periods.  Sub- 
taxa  are  not  additive  to  higher  taxa,  as  reads  are  assigned  according  to  rank-based 
match  percentage  requirements. 
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Figure  14.  Reactor  1  Phylogeny 
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Legend: 

□  Reactor  2  (Early)  ■  Reactor  2  (Acclim)  □  Reactor  2  (O.Oppm) 
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4.4  Analysis 


A  comparison  of  the  changes  in  assigned  reads  between  the  early  and  acclimation 
period  and  the  acclimation  and  late/malathion  addition  period  was  conducted  using 
the  ratio: 

Reads  Assigned^  —  Reads  Assigned^ 

AverageReadsi  —  AverageReads2 

The  genera  most  affected  during  the  transition  from  the  early  period  to  the  acclima¬ 
tion  period  are  shown  in  Tables  9  and  10.  The  genera  affected  by  the  addition  of 
malathion  or  the  later  period  of  operation  (Reactor  2)  are  shown  in  Tables  11  and  12. 
Taxonomic  groups  represented  in  these  figures  are  the  Cytophaga-Flavobacterium- 
Baeteroides  (CFB  group),  Chloroflexi  (green  non-sulfur  bacteria  or  GNSB  group), 
Actinobacteria  (high  GC  gram  positive),  Firmicutes  (low  GC  gram  positive),  Verru- 
comicrobia,  Nitrospirae,  and  Alpha-,  Beta-,  and  Deltaproteobacteria  classes. 

The  first  group  examined  is  the  Proteobacteria  phylum,  which  represented  the 
highest  percentage  group  in  all  three  reactors.  Meganema  is  known  to  cause  bulking 
in  activated  sludge  systems  and  grows  quickly  with  a  high  substrate  affinity  in  high 
oxygen  or  nitrate  environments  (51).  Along  with  Nordella,  Meganema  is  a  member 
of  the  Rhizobiales  order,  which  suggests  a  role  in  nitrogen-fixation.  Bdellovibrio,  a 
predator,  increases  in  high  bacterial  load  environments  and  strongly  affect  biofilm 
composition  (52).  Peredibacter,  another  Bdellovibrionales,  is  also  a  predator  requir¬ 
ing  gram  negative  hosts  in  a  parasitic  cycle  (53).  It  is  possible  that  Bdellovibrio' s 
tendency  to  increase  and  decrease  with  Meganema  in  Tables  9  and  12  indicate  a  preda- 
tor/prey  relationship.  Thauera  is  a  versatile  consumer  of  substrates,  and  can  grow 
aerobically  or  anaerobically  with  nitrate,  which  is  reduced  to  N20  (54).  Thauera 
has  also  shown  a  capability  to  degrade  a  wide  variety  of  aromatic  compounds  (55). 


(6) 
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Table  9.  Increase  after  Acclimation 


Reactor  1  (Acclimation) 

Reactor  2  (Acclimation) 

Reactor  3  (Acclimation) 

Meganema  (causes  bulking) 

Bdellovibrio  (predator) 

Algoriphagus  (NH4  degrader) 

Lewinella  (biopolymer  degrader) 

Meganema  (causes  bulking) 

Leucobacter  (pesticide  degrader) 

Thauera  (aromatic  degrader) 

Lactococcus  (pesticide  degrader) 

Meganema  (causes  bulking) 

Algoriphagus  (NH4  degrader) 

Table  10.  Decrease  after  Acclimation 


Reactor  1  (Acclimation) 

Reactor  2  (Acclimation) 

Reactor  3  (Acclimation) 

Crocinitomix  (possible  OP  degrader) 

Caldilinea 

Peredibacter  (predator) 

Prosthecobacter 

Amaricoccus 

Lactococcus  (pesticide  degrader) 

Mycobacterium  (OPH-like  gene) 

Nitrosomonas  (NH4  to  NO2) 

Nitrospira  (NO2  to  NO3) 

Table  11.  Increase  after  Malathion 


Reactor  1  (O.lppm) 

Reactor  2  (O.Oppm) 

Reactor  3  (3.0ppm) 

Bryobacter  (acid-tolerant) 

Pontibacter  (biopolymer  degrader) 

Cryomorpha  (possible  OP  degrader) 

Caldilinea 

Nordella  (nitrogen  fixation) 

Dyadobacter  (biopolymer  degrader) 

Runella  (biopolymer  degrader) 

Amaricoccus 

Roseomonas 

Caldilinea 

Dysgonomonas  (not  NO3  reducing) 

Nitrospira  (NO2  to  NO3) 

Table  12.  Decrease  after  Malathion 


Reactor  1  (O.lppm) 

Reactor  2  (O.Oppm) 

Reactor  3  (3.0ppm) 

Meganema  (causes  bulking) 

Bdellovibrio  (predator) 

Opitutus  (acid/anoxic) 

Lewinella  (biopolymer  degrader) 

Algoriphagus  (NH4  degrader) 

Nitrospira  (NO2  to  NO3) 

Lewinella  (biopolymer  degrader) 

Meganema  (causes  bulking) 

Leucobacter  (pesticide  degrader) 

Thauera  (aromatic  degrader) 

Lactococcus  (pesticide  degrader) 

Niabella  (biopolymer  degrader) 

Meganema  (causes  bulking) 

Leucobacter  (pesticide  degrader) 

Lactococcus  (pesticide  degrader) 

Leuconostoc  (pesticide  degrader) 

Bdellovibrio  (predator) 
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Americoccus  is  an  aerobic  chemoheterotroph  belonging  to  the  family  Rhodobacter- 
aceae  (56).  Roseomonas  is  an  oxidative  bacteria  associated  with  pathogenicity  and 
infections  (57). 

Algoriphagus,  Dyadobacter,  Pontibacter,  Runella,  Crocinitomix,  Niabella,  Cry- 
omorpha ,  and  Lewinella  are  CFB  group  bacteria,  the  second  largest  representation 
in  the  reactors.  CFB  are  chemoorganotrophic,  generally  aerobic,  and  proficient  in 
degrading  biopolymers  (58).  Algoriphagus  is  a  member  of  the  order  Cytophagales, 
known  for  growth  on  peptone  and  organic  nitrogen  compounds  with  some  groups 
capable  of  degrading  inorganic  nitrogen  compounds  such  as  NH4  (52).  This  may 
have  been  a  factor  in  Reactor  3  where  excess  NH4-N  during  the  acclimation  period 
contributed  to  the  growth  of  Algoriphagus.  Dyadobacter ,  Pontibacter ,  and  Runella 
are  also  members  of  Cytophagales.  Lewinella  is  chemoorganotroph  in  the  Sphingob- 
acteriales  order  and  is  strictly  aerobic  (59).  Crocinitomix  and  Cryomorpha  belong 
to  the  Flavobacteriales  order  (60).  Flavobacterium  sp.  was  the  first  microorganism 
identified  that  could  degrade  OPs  (12).  Cryomorpha ,  in  Table  11,  appears  to  have 
responded  positively  after  the  addition  of  malathion  at  0.1 

Some  of  the  groups  were  represented  by  fermentative  and  acid  tolerant  genera. 
Dysgonomonas,  another  CFB  group  bacteria,  is  fermentative,  facultatively  anaero¬ 
bic,  and  does  not  reduce  nitrate  (61).  Along  with  Nitrospira,  it  showed  an  increase 
in  assigned  reads  after  adding  malathion.  Other  fermentative  bacteria  were  nega¬ 
tively  affected  by  the  addition  of  malathion.  Leucobacter,  along  with  Lactococcus 
and  Leuconostoc,  are  lactic  acid  bacteria  (LAB)  which  are  fermentative  and  prefer 
acidic  environments.  Leucobacter  belongs  to  the  phylum  Actinobaeteria  and  Lac¬ 
tococcus  are  members  of  Firmicutes.  Reactor  1,  which  had  operated  for  over  three 
months,  experienced  a  drop  in  Lactococcus  reads  after  the  transition,  suggesting  that 
a  more  stabilized,  aerobic  nitrification  environment  was  less  suitable  for  fermentative 
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bacteria.  LAB  have  been  shown  to  enhance  the  degradation  of  pesticides  including 
malathion  (62),  but  received  fewer  assigned  reads  during  malathion  addition  in  Re¬ 
actor  3.  This  may  be  the  result  of  a  period  of  adaptation  towards  the  substrate,  or 
could  have  been  affected  by  other  processes.  Leucobacter  and  Lactococcus  were  also 
assigned  fewer  reads  in  Reactor  2,  where  no  malathion  was  added.  Another  genus, 
Opitutus ,  a  Verrucomicrobia  member,  is  also  fermentative  and  believed  to  widespread 
in  anoxic  terrestrial  environments  (52).  The  addition  of  a  foreign  substrate  and  the 
aerobic  conditions  in  Reactor  1  may  have  been  the  cause  of  Opitutus’’ s  decline  after 
malathion.  Finally,  Bryobacter  is  slow-growing,  chemoorganotroph  from  the  phy¬ 
lum  Acidobacteria  that  is  strictly  aerobic  and  also  acid-tolerant  (63).  Bryobacter’ s 
assigned  reads  were  higher  in  Reactor  1  after  the  addition  of  malathion. 

Another  member  of  Verrucomicrobia  is  Prosthecobacter,  an  aerobic  genus  that 
has  been  found  in  ponds,  lakes,  and  sewage  (52).  Caldilinea  is  a  relatively  recent 
discovery,  but  has  been  described  as  a  GNSB-group,  facultatively  aerobic  bacterium 
that  grows  chemoorganotrophically  (64).  Caldilinea  was  assigned  more  reads  in  both 
Reactor  1,  which  had  0.1  ^  malathion,  and  Reactor  2  without  malathion.  Mycobac¬ 
terium, ,  from  the  Actinobacteria  phylum,  have  been  associated  with  tuberculosis  and 
leprosy  (65).  This  genus  also  possesses  the  OPH-like  gene,  that  produces  the  nu¬ 
cleotide  cyclase  phosphodiesterase  Rv0805  (66).  ft  is  not  clear  whether  this  played  a 
role  after  malathion  addition,  based  on  the  number  of  assigned  reads. 

The  remaining  two  genera,  Nitrosomonas  and  Nitrospira,  are  nitrifying  bacteria. 
Nitrosomonas  decline  in  Reactor  1  after  acclimation  is  not  reflected  in  the  nitrifica¬ 
tion  data.  However,  the  total  biomass  had  dropped  to  74%  of  the  higher  mixed  liquor 
TSS  value,  possibly  as  a  result  of  seeding  Reactors  2  and  3.  In  Reactor  2,  during  ac¬ 
climation,  NH4-N  was  low  in  both  the  effluent  and  the  reactor,  suggesting  that  AOB 
were  active.  In  Reactor  3,  effluent  NH4-N  was  present,  suggesting  a  low  population 
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of  AOB.  It  is  likely,  therefore,  that  AOB  were  distributed  primarily  to  Reactor  2.  Ef¬ 
fluent  nitrite  concentrations  were  lower  and  nitrate  concentrations  normal  in  Reactor 
3  during  the  acclimation  period,  corresponding  with  a  decrease  in  assigned  reads  to 
Nitrospira.  This  suggests  that  the  slow-growing,  low  population  of  AOB  was  unable 
to  sustain  the  initial  population  of  NOB,  causing  a  drop  during  acclimation.  The  de¬ 
crease  in  Nitrospira  reads  in  Reactor  1  is  associated  with  the  short-term  rise  in  nitrite 
concentration  six  days  after  the  addition  of  malathion.  However,  in  Reactor  3  where 
the  concentration  of  malathion  was  3  no  similar  decline  in  reads  was  identified  for 
Nitrospira.  Nitrospira  reads  were  higher  for  Reactor  3  after  the  addition  of  malathion. 
Janeczko  (2014)  suggested  that  partial  heterotrophic  decline  may  have  enabled  other 
unaffected  groups  to  increase  and  along  with  those  able  to  metabolize  malathion  (1). 
If  Meganema,  a  species  known  for  bulking,  and  Bdellovibrio,  a  predator,  declined,  as 
shown  in  Table  12,  this  may  have  been  the  case. 
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V.  Conclusions 


This  research  has  provided  insight  into  the  interactions  within  bacterial  commu¬ 
nities  exposed  to  malathion,  a  commercial  pesticide  and  organophosphate.  Diver¬ 
sity  was  highest  in  the  reactor  seeded  directly  from  the  wastewater  treatment  plant 
with  an  initial  Simpon’s  Reciprocal  Index  value  of  15.3.  Temporal  decline  in  di¬ 
versity  was  observed  in  all  three  reactors.  Ammonia-oxidizing  bacteria  exposed  to 
malathion  for  30  days  at  0.1  y^-  and  3.0  '-y-  were  not  significantly  affected.  Nitrite- 
oxidizing  Nitrospira  was  associated  with  increased  reads  during  exposure  to  0.1  y1 
and  decreased  reads  during  exposure  to  3.0  y1.  The  genus  Cryomorpha,  a  mem¬ 
ber  of  known  organophosphate-degrading  order  Flavobacteriales,  was  associated  with 
increased  reads  after  30  days  exposure  to  malathion  at  0.1  y1. 
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